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ABSTRACT 


The  classic  beamforming  algorithm  involves  the  process  of  "phasing- 
to-a-plane."  This  report  presents  a method  of  beamforming  which 
does  not  require  the  phasing-to-a-plane  step,  and  therefore 
can  be  implemented  without  the  use  of  time  delay  or  phase  shift 
networks . 

The  study  was  restricted  to  cylindrical  or  circular  arrays  of 
equally  spaced  sensor  elements,  although  the  method  appears  to 
be  equally  applicable  for  other  geometries.  Since  no  general 
closed  form  solution  for  the  required  shading  coefficients  is 
available,  a set  of  shading  coefficient  selection  algorithms  are 
developed  and  presented  along  with  the  resulting  beam  pattern 
performance.  Studies  of  the  bandwidth,  elevation  beam  patterns, 
and  the  effects  of  construction  errors,  as  well  as  methods  of 
implementing  the  beamformer,  are  reported. 

The  results  of  the  study  and  the  experimental  program  indicate 
that  substantially  equal  performance  can  be  achieved  with  a 
physically  smaller  array.  In  addition,  the  correlation  method 
is  easier  to  implement.  As  an  example,  the  Two  Degree  Array 
(TDA)  was  compared  to  the  Cylindrical  Transducer  (CT)  for  an 
experimented,  sonar.  The  TDA  was  only  65#  the  size  of  the  CT.  The 
TDA  contained  the  same  number  of  elements  spaced  in  the  azimuthal 
direction,  although  it  contained  l/9  as  many  elements  in  the  elevation 
direction:  the  TDA  and  CT  operated  at  the  same  frequency  and 

demonstrated  essentially  equal  performance  in  both  the  azimuthal  and 
elevation  planes.  Additional  examples  of  size  and  cost  advantage 
are  given  in  Appendix  C. 

The  chapter  on  array  rise  time  and  the  use  of  coded  pulse 
transmission  to  accomplish  pulse  compression  on  the  array  is 
equally  applicable  to  phase  shift  beamformers . 

The  assorted  effects  of  changing  the  sector  of  active  elements 
on  the  cylinder  are  presented  in  Appendix  A.  Appendix  B reports 
the  results  of  comparing  the  algorithmic  approach  to  an  equiva- 
lent design  by  using  linear  programming  as  an  optimization 
technique . 

Appendix  C is  a design  analysis  of  the  reduction  in  size  and  cost 
savings  to  be  expected  by  using  the  correlation  beamformer  method 
as  applied  to  three  existing  sonar  systems.  Appendix  D is  a repro- 
duction of  a signal- to-nolse  ratio  study  comparing  the  correlation 
beamformer  with  a phase  shift  beamformer. 
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W - the  subtended  width  of  the  p - distance  from  array  center  to 

array  as  seen  from  a point  in  field  point  in  wavelengths , 

the  farfield  on  the  DMRA  p = £ . 

A. 

(normalized  to  the  design  ^o’^o’^o’  etc*  ” des*6n  v®1-11®13 

wavelength);  alternately  used  the  variables, 

to  designate  the  line  proper,  t - the  phase  term  contained  in  A^; 
i.e.,  "the  W line".  the  exact  phase  of  the  signal 

on  the  array  at  the  selected 
time  of  maximum  response; 
introduced  p.  11. 
o>  - the  angular  frequency  of  a 
signal  in  radians  per  second. 
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I.  INTRODUCTION 


In  January  1973,  a method  for  forming  beams  using  elements  or  staves 
on  a cylindrical  array  without  the  use  of  time  delays  or  phase  shift 
networks  was  discovered.  Relatively  large  arrays  were  under  consideration, 
that  is,  cylinders  whose  radii  in  wavelengths  were  greater  than  about  10. 

A quick  survey  of  the  literature  revealed  no  additional  information  on  the 
method,  and  further  computer  simulation  uncovered  no  apparent  problems  or 
serious  disadvantages  that  might  result  from  using  the  processing  scheme. 

The  Naval  Coastal  Systems  Laboratory  was  asked  to  support  ARL 
during  a two  year  study  of  the  beamforming  method.  This  report  is  the 
major  result  of  that  study. 

During  the  first  year,  a computer  model,  which  was  used  throughout 
the  remainder  of  the  study,  was  developed,  and  two  experiments  were 
performed  to  validate  the  model.  A partial  understanding  of  the  role  of 
the  coefficients  and  their  effect  upon  the  beam  pattern  performance  was 
gained,  and  the  problem  of  rise  time  was  delineated.  In  addition,  the 
realization  that  a beam  could  be  scanned  at  high  speed  began  to  take 
shape . 

During  the  second  year,  a study  of  high  speed  scanning  and  its 
implementation  was  begun,  and  other  methods  of  implementing  preformed 
beams  were  studied.  The  pulse  compression  problem  was  studied  and  a 
pulse  code  was  developed  which  achieves  pulse  compression.  Larger 
arrays  were  studied  for  possible  systems  applications  and  for  preliminary 
design  of  hardware  to  be  used  in  testing  the  pulse  compression  codes.  A 
more  complete  understanding  of  the  role  of  the  coefficients  was  obtained. 
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The  purpose  of  this  report  is  two- fold:  (l)  to  document  the 

essence  of  the  correlation  beamforming  method  and  to  summarize  what 
was  learned  during  the  course  of  the  study,  and  (2)  to  serve  as  a 
design  guide  for  those  who  would  use  the  method. 

Chapters  III  and  IV  of  this  report  describe  the  fundamental 
concepts  of  the  correlation  beamformer  method,  and  describe 
"rules-of- thumb"  that  can  be  used  to  design  or  configure  a 
correlation  array.  One  important  appendix  is  C,  which  discusses  the 
cost  and  size  advantage  to  be  obtained  by  the  use  of  the  method. 


II.  BACKGROUND 


Beam  formation  is  a well-known  process  by  which  the  directionality 
of  response  of  an  array  of  acoustic  transducers  is  controlled.  In  the 
general  case,  the  individual  transducer  elements  in  an  array  may  be 
arranged  with  any  spatial  distribution,  and  may  have  their  individual 
response  maxima  oriented  in  any  direction. 

Beam  formation  can  be  accomplished  in  principle  in  the  general  case; 
however,  only  in  the  simplest  cases  can  a closed-form  solution  for  the 
shading  coefficients  be  obtained.  If,  for  example,  the  elements  are 
equally  spaced  along  a line  and  either  have  no  response  maxima  (omnidirec- 
tional), or  are  all  pointed  in  the  same  direction,  then  one  can  solve  for 
the  shading  coefficients  so  that  the  beam  pattern  will  have  specific 
properties.  The  literature  on  beam  control  is  extensive  in  this  case. 

Another  relatively  simple  case  consists  of  a set  of  omnidirectional 
elements  equally  spaced  on  a circle;  there  is  no  baffle  so  that  all 
elements  can  be  seen  from  all  directions.  Some  literature  for  this  case, 
which  contains  solutions  for  the  shading  coefficients,  is  available,  such 
as  Ref.  1. 

For  arrays  more  complex  than  the  simple  ones  as,  for  example,  a 
cylindrical  baffle  containing  finite  sized  elements  pointing  in  a radial 
direction,  no  closed-form  solutions  for  the  shading  coefficients  have  been 
obtained.  In  such  a case,  the  classical  method  of  "phasing  to  a plane" 
is  used. 

Conventional  beeunforming  or  spatial  processing  is  accomplished  through 
a linear  weighted  combination  of  the  outputs  of  the  elements.  That  is, 
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the  beamformer  output  is 


N 

E*(cp,t)  = £ A*  e*((p,t)  , (2-1) 

n=-N 

where  e*(<p,t)  1b  the  complex  voltage  developed  by  the  nth  element  when 
in6onified  by  the  signal  S(t)  arriving  from  the  angle  qp  and  at  time  t, 
and  A*  is  the  complex  weighting  coefficient  which  accomplishes  beam 
formation  (both  phasing  and  shading) . 

The  problem  is  to  select  the  A*  in  such  a way  that  the  envelope  of 
E*(qp,t)  has  one  maximum  value  at  angle  cpQ  and  tQ,  and  has  controlled  low 
levels  at  other  angles  q#cpo,  and  other  times  t^tQ. 

The  classic  algorithm  involves  the  selection  of  A*=A^exp( i0^)  in  two 
parts.  First,  0r  is  selected  to  "phase-the-element-to-a-plane"  or  delay 
its  output  by  some  angle  0^,  so  that  the  collection  of  M elements  appear 
to  lie  in  a line  or  on  a plane.  For  a time  delay  beamformer,  0R  may 
involve  several  revolutions.  For  a phase  shift  beamformer,  0R  is 
restricted  to  ±it. 

The  second  part  of  the  classic  algorithm  is  the  selection  of  the 
magnitude  of  the  shading  coefficients,  the  Aq.  This  operation  proceeds 
according  to  methods  contained  in  the  extensive  literature  concerning 
line  array  shading  coefficients  (for  example,  see  Kefs.  2 through  6), 
and  then  computer  tests  are  run  to  obtain  the  final  result. 

In  both  time  delay  and  phase  shift  beamforming,  the  angle  0^  is 
important,  and  during  implementation  requires  the  use  of  phase  shifters 
or  delay  lines.  In  sonar  arrays  of  high  resolution  covering  a wide  sector, 
such  devices  preclude  scanning  a single  beam  because  of  the  restricted 
data  rate,  and  can  become  a major  expense  item  when  many  preformed  beams 
are  implemented  in  order  to  meet  the  data  rate  requirements. 
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During  late  1972,  while  performing  work  on  the  suppression  of 
grating  lobes  in  the  case  of  widely  spaced  elements  on  cylindrical  arrays, 
it  was  discovered  that  beams  could  be  formed  without  the  "phasing-to-a- 
plane"  step.  If  useful  beams  could  be  formed  and  controlled  without  the 
phase  shifters  or  delay  lines,  then  arrays  of  smaller  size  at  less  expense 
could  be  constructed.  Therefore,  a proposal  to  study  the  "method  of 
beamforming  without  the  use  of  time  delay  or  phase  shift  networks"  was 
submitted  and  accepted.  This  report  is  the  result  of  that  study. 

The  cylindrical  array  was  chosen  for  study  because  of  its  azimuthal 
symmetry;  many  beams  can  be  formed  by  applying  the  same  set  of  shading 
coefficients  to  a different  set  of  elements. 

The  method  of  "forming  beams  without  the  use  of  time  delay  or  phase 
shift  networks"  or  "forming  beams  by  the  use  of  real  shading  coefficients" 
has  finally  come  to  be  called  a "correlation  beamformer, " primarily 
because  of  the  way  the  coefficients  are  selected.  However,  the  new  method 
does  not  simplify  the  mathematics  of  beam  formation;  therefore,  as  before, 
no  closed- form  solution  is  available. 

The  solution  to  the  problem  of  finding  shading  coefficients  for  the 
correlation  beamformer  on  a cylindrical  array  has  been  to  develop  a set 
of  rules  or  algorithms  to  determine  the  coefficients;  then  a computer 
implemented  mathematical  model  is  exercised  to  evaluate  the  predicted 
results.  By  this  process,  the  heuristically  derived  coefficient  selection 
algorithms  have  been  proved,  and  array  scaling  rules  have  been  determined. 

Chapters  III  and  IV  discuss  the  array  performance  in  relation  to  the 
coefficient  selection  algorithms.  The  array  fill  time  problem  and  two 
proposed  solutions  are  discussed  in  Chapter  V.  Chapters  VI,  VII,  and  VIII 
of  the  report  include  methods  of  implementing  either  scanned  or  preformed 
beams,  some  technological  fallout  of  the  study,  and  experimental  verifi- 
cation, respectively.  Chapter  IX  gives  the  study  conclusions  and 
recommendations  for  further  work. 
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Although  this  study  was  aimed  at  the  correlation  beamformer  for 
cylindrical  arrays,  the  correlation  method  is  equally  applicable  to 
other  array  geometries  (for  example,  lines  and  planes  or  volume  arrays). 
The  coefficient  selection  algorithms  are  unchanged  for  other  geometries. 
Also,  some  of  the  performance  results  determined  in  this  study  apply 
equally  well  to  phase  shift  or  time  delay  beamforming  methods. 
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III.  THE  DESIGN  AND  PERFORMANCE  OF  THE  REFERENCE 
CYLINDRICAL  ARRAY- -THE  SIMPLE  CASE 


A.  The  Role  of  the  Shading  Coefficients 

In  this  chapter,  the  study  of  the  design  and  performance  of  the 
correlating  beamformer  will  be  restricted  to  the  simplest  possible  case: 
point  elements,  equally  spaced  one-half  wavelength  apart  on  the  front 
half  of  a cylindrical  baffle  of  radius  R wavelengths.  The  coefficient 
selection  algorithm  will  be  described,  and  the  beam  performance  and 
scaling  rules  developed. 

Figure  3-1  shows  the  assumed  geometry.  For  the  present,  the  test 
probe  will  be  assumed  to  be  in  the  farfield  and  will  emit  a simple  CW 
signal  sin(cnt).  The  baffle  will  be  assumed  to  be  opaque  but  otherwise 
does  not  affect  the  element's  pattern;  that  is,  an  element's  response  is 
constant  and  independent  of  angle  if  it  is  in  the  line  of  sight  of  the 
test  probe;  the  output  is  zero  otherwise. 

The  beamformint  operation  will  be 

N 

E*(<P,t)  = E An  , (3-1) 

n=-N 


where  the  * denotes  a complex  number  and  e*  iB  the  output  voltage  from  the 
nth  element  at  some  time  t when  the  test  probe  is  at  angle  <p. 


WM* 


where 


A-^n)  is  the  signal  replica  factor, 

Ag(n)  is  the  sensor  density  corrector,  and 
Aj(n)  is  the  aperture  window  function. 

1.  The  Signal  Replica  Factor 

Consider  a cylinder  whose  radius  is  15  wavelengths.  Further, 
assume  the  n=0  element  is  in  the  direction  of  the  major  response  axis 
(MRA)  (the  X axis  in  Fig.  3-l).  If  the  test  probe  is  at  qfeO  and  the  array 
sensitivity  is  measured,  the  result  will  be  one  apparent  element.  As 
additional  elements  are  added  in  pairs,  n=±l,±2,  etc.,  the  array  sensitivity 
will  grow  until  the  additional  elements  lie  outside  the  first  Fresnel  zone 
on  the  cylinder.  As  elements  are  added  into  the  second  Fresnel  zone,  the 
array  sensitivity  will  decrease,  making  the  apparent  number  of  elements 
on  the  array  decrease.  The  dotted  curve  in  Fig.  3-2  illustrates  this 
behavior. 

If  phase  shift  or  time  delay  beamforming  were  used,  the  elements 
would  be  moved  electrically  so  that  they  all  would  add  coherently,  and  the 
dashed  curve  would  result.  One  can  see  from  the  dotted  curve  that  the 
addition  of  the  pair  of  elements  marked  15  will  cause  the  overall  sensi- 
tivity of  the  array  to  decrease  because  these  added  elements  exist  in  the 
second  Fresnel  zone.  If  one  were  to  reverse  their  signs  before  adding, 
the  solid  curve  would  result.  That  is,  the  solid  curve  is  the  result  of 
weighting  by  +1  those  elements  in  the  odd  numbered  zones,  and  weighting 
by  -1  those  elements  in  the  even  numbered  zones. 

The  choice  of  A^(n)=±l  does  result  in  a narrow  major  lobe  with 
good  sensitivity,  but  with  unacceptable  side  lobe  performance.  The  concept 
of  a replica  correlator  suggests  that  the  factor  A^n)  be  selected  to 
more  accurately  replicate  the  signal  at  an  element  location. 
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Figure  3-3  shows  a long  CW  signal  insonifying  an  array  of 
elements.  In  order  to  define  a signal  replica,  we  must  know  the  signal 
type,  the  proposed  direction,  and  the  proposed  times  at  which  a peak 
positive  beamformer  output  is  desired.  For  example,  assume  a simple  sine 
wave  arriving  from  p=0,  and  require  a peak  output  at  the  instant  when  the 
signal  is  a cosine  vave  on  the  array  as  shown  in  Fig.  3-5 • The  nth 
element  at  angle  a from  the  MRA  or  signal  matching  direction  is  shown. 

The  nth  element  is  R-R  cos  0^  wavelengths  behind  the  n=0  element  in  the 
direction  of  signal  propagation.  So 

2xR(l  - cos  a^)  (3-5) 

represents  the  distance  in  radians  beyond  the  zeroth  element.  Thus,  the 
impressed  instantaneous  signal  at  element  n is 

cos[2rtR(l  - cos  Qin)]  (5-1*-) 

provided 

laJSn/2  . (3-5) 

If  the  signal  replica  is  modified  with  a phase  term  \|r,  such  that 

Ai(an)  = cos[2xR(l  - cos  aR)  - \|r]  , (3-6) 

then  it  can  be  made  to  represent  any  desired  signal  phase  shift  before 
matching.  For  the  simple  case  under  consideration,  however,  \(r= 0,  as  Bhcwn 
in  Fig.  3~3>  and  the  signal  replica  factor  is 

A^(a  ) = cosl2jtR(l  - cos  an)]  . (3-7) 

Chapter  IV  contains  a section  devoted  to  the  more  general  case  of 
arbitrary  \|r. 
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2 . Sensor  Density  Corrector 


If  the  coefficient  factor  A^(n)  only  is  used,  and  a signal  is 
applied  to  the  array,  one  would  in  general  find  a varying  sensitivity 
over  the  aperture  W shown  in  Fig.  3-1.  Of  course,  the  varying  sensitivity 
over  the  W aperture,  called  the  aperture  window  function,  is  usually  used 
to  control  side  lobe  levels  in  line  or  plane  arrays.  This  window  function 
must  have  special  characteristics  if  it  is  to  be  useful  in  controlling 
side  lobe  levels.  Our  approach  to  this  problem  is  in  two  parts:  (l)  a 

function  A^n)  which  will  make  the  A^(n)  * A^(n)  aperture  window  function 
look  uniform,  and  (2)  a function  A^(n)  which  independently  applies  the 
required  shading  for  side  lobe  control. 

The  method  of  selecting  Ag(n),  the  sensor  density  corrector,  is 
to  first  evaluate  the  apparent  sensitivity  over  the  W aperture  when  A^n) 
only  is  applied.  One  of  the  main  contributors  to  the  aperture  function  is 
the  apparent  density  of  sensors  when  they  are  projected  to  the  W line. 

Other  contributors  can  be  the  envelope  of  A^  or  the  envelope  of  the  signal. 
We  are  not  interested  in  the  cyclic  nature  of  A^  or  the  carrier  of  the 
signal  over  the  aperture,  only  the  average  sensitivity  per  unit  length  on 
W;  therefore,  only  the  envelopes  are  used. 

As  an  example,  consider  some  unit  width  on  W near  the  center  of 
the  array.  The  signal  strength  there  is  unity.  Also,  the  envelope  of 
A1(a)  is  unity.  Further,  the  projected  sensor  density  is  about  the  same 
as  the  density  on  the  cylinder,  which  is  two  elements  per  wavelength. 

Now,  consider  the  situation  near  the  end  of  W.  Again,  the  signal  and  A^ 
envelopes  are  unity,  but  the  projected  sensor  density  near  the  edge  of  the 
array  is  very  large.  The  aperture  window  function  is  therefore  larger 
toward  the  edges  of  the  aperture,  and  A^a)  should  be  chosen  to  correct 
this  effect. 

In  Fig.  3-1,  let  the  distance  along  the  W line  be  Y,  with  Y=0 

at  the  center  of  the  W line.  An  element  at  a on  the  cylinder  is 

J 

15 


a 


projected  onto  the  W line 


at  YJf 


where 
= R sin  a 

J 


If  the  element  dens' ty  on  the  cylinder  is 

q _ number  of  elements 
da  unit  angle 

then  the  projected  density  on  W is 


(3-8) 


(3-9) 


.S.  = _! i 

dY  da  R cos  a 


(3-10) 


Since  q/Rda  is  a constant,  the  functional  dependence  of  element  density  on 
the  W line  is  l/cos  a.  The  obvious  choice  for  the  density  corrector,  then, 
is  A2(a)=cos  a,  since  the  aperture  function  for  the  product  A^(a)  * A^(a) 
will  be 


■£7  = "5 — - — — cos  a = constant  . (3-H) 

dY  da  R cos  a 

If  the  signal  envelope  had  varied  over  the  array,  then  A^(a)  would  have 
varied,  and  these  effects  would  have  been  incorporated  in  the  calculation. 

In  fact,  in  Chapter  V,  a triangular  envelope  signal  will  be  used  and 
Ag(a)  calculated  for  that  case. 

3.  The  Aperture  Window  Function 

This  part  of  the  coefficient  selection  algorithm  is  based  on  the 
approximation  that,  for  beam  pattern  properties  near  the  major  lobe 
(beamwidth,  first  side  lobe  level,  sensitivity),  the  array  aperture 
function  appears  to  be  approximately  independent  of  field  angle  cp.  It  is 
further  assumed  that,  for  angles  near  the  MRA,  the  Fourier  transform 
relationship  exists  between  the  array  illumination  function  and  the  farfield 
directional  response,  as  it  does  for  a line  array. 


Now,  we  have  a signal  replica  and  an  A 2 selected  so  that  the 
window  function  is  uniform.  Line  array  theory  would  predict  a first  side 
lobe  level  of  -13  dB//maJor  lobe  level.  If  we  contrive  by  the  faccor 
A^(a)  to  make  the  aperture  function  on  W some  function,  say  F(Y),  then  by 
analogy  with  line  arrays,  the  farfield  directivity  pattern  should  be 
described  by  the  Fourier  transform  of  F(Y). 


For  example,  let  A^(a)  be  arbitrarily  selected  as 
A^(a)  = cos  a 


(3-12) 


The  projection  of  this  function  on  the  W line  can  be  found  by  noting  that 

Y = R sin  a , (3-13) 


or 


a = arcsin  ^ 


(3-1*0 


If 

then 

or 


A^(a)  = cos  a , 


A(Y)  = cos(arcsin  ^) 


(3-15) 

(3-16) 

(3-17) 


The  normalized  modulus  of  the  Fourier  transform  of  the  A(Y)  can  be  shown 
to  be 


D(tp) 


2 Jx(x) 

x 


(3-18) 


where  x = 2xR  sin  (p  and  is  a Bessel  function  of  the  first  kind  of 
order  one.  (See  Ref.  7.) 
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The  A^(a)  algorithm  states  that  the  function  D(cp)  above  should 
be  descriptive  of  the  beamwidth  and  first  side  lobe  level  if  the  A^ 
conjecture  is  nearly  true.  The  function  D(sp)  is  well  known  and  is  the 
same  as  the  response  of  a broadside  circular  piston  of  radius  R wave- 
lengths. The  half-power  beamwidth  of  the  function  is 

BW  = 2 arcsin  (radians)  , (3-19) 

and  the  first  side  lobe  is  suppressed  1J.8  dB  below  the  level  of  the 
major  lobe. 

A simple  development  of  the  Fourier  transform  relationship  for 
line  and  plane  arrays  is  given  in  Ref.  7 and  an  excellent  treatment  of 
Fourier  integral  methods  in  beam  pattern  analysis  is  given  in  Ref.  4-. 

The  final  shading  coefficients  A(a^) , represented  by  the  product 
of  its  factors  as  shown  in  Eq.  (3-2)  for  the  simple  reference  array,  are 

O 

A(a)  = cos  a cos[2nR(l  - cos  a)]  . (3-20) 

This  continuous  function  of  a is  sampled  at  the  locations  of  each  element 

a to  arrive  at  the  list  of  discrete  shading  values  to  be  used  in  the 
n 

beamformer. 

Figure  3-^  shows  plots  of  Eq.  (5-20)  for  two  small  values  of  R. 
The  figure  should  help  visualize  the  shape  of  the  shading  coefficients  and 
their  dependence  on  the  radius  parameter  R. 

B.  The  Performance  of  the  Array 

The  reference  array  is  an  array  of  point  sensors  spaced  \/2  apart, 
covering  the  front  half  of  the  surface  of  a cylindrical,  opaque  baffle  of 
radius  R wavelengths.  The  test  signal  is  an  infinitely  long  CW  signal, 
sin(<nt),  emitted  from  a probe  in  the  very  farfield  at  angle  <p.  The 
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selected  shading  coefficients  sire,  as  developed  in  the  previous  sections. 


A(an)  = cos2an  cos[2nR(l  - cos  otn) ] . (5-21) 

The  computer  model  was  exercised  several  times  for  this  case,  using 
the  value  of  R as  the  prime  variable,  in  order  to  study  the  effects  of 
scaling. 


1.  The  Azimuthal  Beam  Pattern 

Figure  5-5  shows  the  azimuthal  beam  pattern  for  R-12.5  and 
N=38(2N+l=M=total  number  of  elements).  The  most  striking  observation  is 
that  the  side  lobe  region  does  not  show  the  character  displayed  in  patterns 
of  phase  shift  or  time  delay  beamformers.  There  is,  instead,  an  almost 
constant  level  or  "pedestal"  of  sensitivity  nearly  independent  of  angle. 

We  have  defined  the  "average  side  lobe  level"  (SLL)  as  an  rms  average 
sensitivity  of  the  side  lobe  region  from  the  first  null  near  the  major 
lobe  out  to  some  final  angle,  usually  60°  in  pattern  of  the  type  Bhown  in 
Fig*  5-5*  Although  the  term  SLL  implies  absolute  level,  we  actually 
present  "average  side  lobe  suppression,"  an  average  suppression  below  the 
major  lobe,  in  decibels.  Thus,  the  SLL  of  Fig.  3-5  is  18.9  dB. 

Several  different  schemes  for  selecting  shading  coefficients 
were  tried  in  an  attempt  to  change  the  SLL  of  Fig.  3-5*  After  a rather 
extensive  study,  in  which  only  the  variance  or  ripple  of  the  side  lobe 
structure  could  be  changed,  it  was  concluded  that  this  SLL  pedeBtal  appears 
to  be  essentially  unaffected  by  any  design  parameter  other  than  array 
size  R, 

This  constant  side  lobe  level  or  pedestal  behavior  has  been 
noted  by  Steinberg  (Ref.  8)  for  line  arrays  of  randomly  spaced  elements. 

His  explanation  for  line  arrays  is  simply  that,  on-axis,  all  elements  see 
the  same  pressure,  and  hence  the  summation  in  the  processor  (adder)  is  a 
coherent  addition.  As  more  elements  are  added  to  the  array,  the  on-axis 
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sensitivity  increases  as  M increases,  M being  the  number  of  elements. 

For  off-axis  pressure  (i.e.,  the  side  lobe  region),  the  elements  see 
signals  that  are  random  in  phase  (because  of  the  random  spacing)  and 
the  processor  effects  an  incoherent  addition,  that  is,  side  lobe  sensi- 
tivity increases  as  ^M.  The  side  lobe  level  with  reference  to  the  MRA 
level,  therefore,  is  given  by 

-20  log  = -10  log  M . (3-22) 

\J  M 

Steinberg  derives  the  mean  power  pattern  of  a randomly  spaced 
line  array  as 

|f(U)|2  = |f0(u)|2(l  - jjj)  + i . (3-23) 

He  explains  that 

"The  first  term  is  the  desired  power  pattern*  slightly 
reduced  in  strength  to  account  for  the  angle  independent 

Tedestal  l/M.  In  the  neighborhood  of  the  main  lobe 
f0|2»l  and  rises  well  above  the  pedestal.  Hence,  the 
shapes  of  the  main  lobe  and  of  those  nearby  side  lobes 
whose  strength  »l/M  are  essentially  unchanged.  Away 
from  the  main  lobe,  where  the  side  lobe  level  of  the 
design  pattern  <l/M,  the  character  of  the  design  pattern 
is  submerged  below  the  random  component  and  disappears." 

Another  point  of  significance  in  Steinberg's  paper  is  the 
conclusion  that  deterministic,  aperiodically  spaced  arrays  and  random, 
aperiodically  spaced  arrays  that  he  considered  showed  statistically 
indistinguishable  behavior  of  the  side  lobes. 

Thus,  the  suggestion  was  strong  that,  although  the  elements  of 
the  cylindrical  array  are  equally  spaced  on  the  cylinder,  they  are  rather 
aperiodically  spaced  in  the  direction  of  an  approaching  wavefront  and 
might  therefore  meet  one  of  Steinberg's  conditions. 

i 

fQ,  we  think,  is  controlled  by  A^(a). 


d 
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Further,  the  beam  patterns  that  result  with  the  cylinder  are 
rather  well  described  by  the  Steinberg  approach.  For  example,  the  Steinberg 
pedestal  should  be  10  log  M below  the  level  of  the  MRA.  Since  the  cylindri- 
cal array  of  radius  F wavelengths  has  a circumference  of  2itR,  and  since 
there  are  elements  spaced  each  \/2  on  one-half  the  array,  then  there  are 
2rtR  elements  in  each  aperture. 


This  conjecture  was  tested  for  arrays  with  R=12  to  R=407  in 
steps  of  approximately  a factor  of  two.  In  every  case, 

SH  = -10  log  2*R  (5-24) 

described  the  resulting  average  side  lobe  suppression  to  within  1 dB. 
Figure  5-6  is  an  example  of  a beam  pattern  plot  with  R=100.  The  pre- 
dicted SLL  is  -28  dB. 


Steinberg's  average  side  lobe  level  was  derived  assuming  equal 
strength  elements,  that  is,  no  shading,  or  A(y)=l.  If  the  elements  were 
shaded,  for  example,  A^*l,  then  the  on-axis  sensitivity  would  become 

M 

E A , (3-25) 

i=l 

not 

M 

E 1 = M > (5-26) 

i=l 


and  the  incoherent  addition  in  the  side  lobe  region  would  become 


not 


(5-27) 


(5-28) 
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In  spite  of  these  differences,  which  have  so  far  remained  unexplained, 
the  prediction 


SLL  = -10  log  2nR  = -10  log  M (3-29) 

appears  to  be  good. 

In  the  Steinberg  quote  on  a previous  page,  the  statement  that 
"those  features  of  the  desired  pattern  whose  strength  is  much  greater 
than  l/M  are  not  changed  by  the  presence  of  the  pedestal"  would  suggest 
that  the  beamwidth  (BW)  of  the  pattern  should  be  determined  by  the  D(cp) 
function,  derived  from  the  transform  of  the  A^(Y). 

Now,  when 


then 


and 


where 


A^(a)  = cos  a 


D(cp) 


2 Jx(x) 
x 


x = 2nR  sin  cp 


Solving  for  cpQ=BW/2  in  the  equation 


leads  to  the  result  that 


BW  = 2 arcsin  (ra<*lans) 


(3-30) 

(3-31) 

(3-32) 

(3-33) 

( 3-3*0 

(3-35) 
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Again,  this  prediction  was  compared  to  the  "measured"  beamwidth 
of  the  computed  pattern  from  R=12  to  R=407,  with  the  result  that  Eq.  (3—35) 
predicts  beamwidth  to  within  2$. 

Another  result  of  the  beamwidth  computations  is  that  the  "beam 
spreading"  is  less  than  20$.  Beam  spreading  is  defined  as  the  percentage 
increase  in  resulting  beamwidth  over  that  beamwidth  which  can  be  expected 
with  a uniformly  illuminated  aperture  of  the  6ame  length  (2R) . It  is  a 
measure  of  beam  spreading  caused  by  the  application  of  some  shading 
function.  For  example,  cosine  squared  shading  on  a line  array  spreads 
the  beamwidth  64$  over  that  attainable  on  the  same  aperture  with  uniform 
shading.  The  20$  figure  obtained  with  the  present  method  indicates  a 
rather  large  utilization  of  available  aperture. 

Given  a shading  function,  in  this  case  A (a)=cos  a,  one  can 
define  an  aperture  constant  K,  which  is  the  product  of  the  number  of 
wavelengths  in  the  aperture  times  the  resulting  half-power  beamwidth  in 
degrees,  provided  the  BW  is  not  too  large.  Thus 

K = 2R  • BW  (degrees  • wavelengths)  . (3-36) 

In  the  present  case  (the  reference  array) 

K = 2R  • ^2  . 2 arcsin  . (3-37) 

When  R>10,  BW(deg)<3°,  and  sin  <P0»CP0  to  within  less  than  0.1$,  Eq.  (3-37) 
can  be  approximated  by 

K = ^ * 59(  wave  length  • degrees)  . (3-38) 

it 

An  additional  observation  concerns  the  behavior  of  the  first 
side  lobe,  usually  the  largest  lobe,  which  is  indicated  by  SLI^.  Recall 
that  since  A^(a)=cos  a,  the  predicted  coherent  beam  pattern  is  D(qp)  as 
given  in  Eqs.  (3-32)  and  (3-33)*  The  suppression  of  the  first  side  lobe 
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of  this  function  Is  17.8  dB,  independent  of  array  size.  However,  the 

incoherent  pattern  or  random  pedestal,  which  does  depend  on  array  size, 

may  add  or  subtract  from  the  predicted  peak  side  lobe  level.  In  Fig.  3-5 

the  random  pedestal  is  at  about  -19  dB.  The  resulting  lobe,  considering 

the  addition  or  subtraction  of  the  two  components  with  unknown  phase, 

could  range  from  a high  of  -12  dB  to  a low  of  -3 6 dB.  So  it  is  not 

surprising  that  the  peak  lobe  in  Fig.  3-5  does  not  measure  -17-8  dB. 

The  same  argument  can  be  applied  to  Fig.  3-6.  The  range  of  values  of 

peak  side  lobe  level  that  should  not  be  surprising  in  Fig.  3-6  are  from 

a high  of  -15.2  dB  to  a low  of  -21.5  dB.  The  main  point  concerning  the 

peak  side  lobe  SLL  is  that  it  is  fruitless  to  choose  an  A, (Y)  which 
P 3 

predicts  a side  lobe  less  than  -10  log  2nR,  since  it  would  Just  disappear 

below  the  random  component. 

2.  The  Elevation  Beam  Pattern 


Figure  3-7  shows  the  elevation  beam  pattern  for  two  sizes  of 
the  reference  case  array.  Several  such  plots  were  studied,  and  all 
behaved  as  shown.  The  half- power  beamwidth  decreases  with  the  inverse 
square  root  of  size,  much  like  the  beamwidth  of  an  end-fire  array.  In 
fact,  the  measured  beamwidth  from  computer  runs  can  be  predicted  from  the 
end- fire  array  directivity  function  if  one  uses  a length  of  0.55  R.  For 
the  reference  case,  the  beamwidth  can  be  calculated  from 

BW(e)  = -^  (deg)  . (3-39) 

This,  of  course,  depends  heavily  upon  the  shading  coefficients  and 
especially  the  sector  of  active  elements,  in  this  case  2a^=l80°. 

The  side  lobe  region  shown  in  Fig.  3-7  is  very  unlike  an  end-fire 
array,  however.  No  clear  explanation  of  the  smooth  behavior  in  Fig.  3-7 
is  available.  The  sensitivity  at  6=90°,  that  is,  with  the  cylindrical 
array  seen  end  on,  is  a measure  of  the  average  value  of  the  coefficients. 
This  is  because,  at  this  angle,  the  signal  insonlfying  the  array  is  in 
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phase  at  each  element.  The  suppression  at  6=90°  has  been  found  empirically 
to  be  roughly  described  by 

-10  log  R yiF  dB//major  lobe  level  . O-^+O) 


3-  Bandwidth — Frequency  Response 

Again,  because  the  array  has  depth  in  the  direction  of  signal 
propagation,  signal  detection  is  expected  to  be  a rather  narrowband 
process.  During  the  matching  operation,  that  is,  during  the  selection 
of  the  A^(an)  coefficients,  R cycles  of  the  signal  must  be  known.  If, 
for  example,  during  operation  of  the  array,  the  center  frequency  was 
changed  so  that  R+l/2  cycles  of  the  signal  occurred  in  the  original  radius 
of  the  array,  then,  when  the  front  edge  of  the  array  was  insonified  by  the 
start  of  a cosine  wave,  the  back  of  the  array  would  be  insonified  by  a 
signal  exactly  out-of-phase  with  what  was  expected. 


In  order  to  test  the  array  bandwidth  behavior,  the  computer 
model  was  run  for  several  operating  frequencies  for  an  R=I2.3  array  chosen 
to  operate  at  100  kHz.  Figure  3-8  is  a summary  of  these  runs,  and  shows 
how  the  various  pattern  parameters  change  with  operating  frequency. 


If  the  bandwidth  is  considered  as  "that  frequency  range  over 
which  the  array  sensitivity  has  deteriorated  no  more  than  3 dB, " then  the 
bandwidth  is  14.5  kHz.  However,  over  this  range,  the  beamwidth  increases 
by  a factor  of  two!  If  the  bandwidth  is  restricted  to  that  range  of 
frequencies  over  which  the  beamwidth  does  not  increase  beyond  25%,  then 
the  operating  bandwidth  becomes  7*3  kHz.  If  the  bandwidth  is  defined  as 
"the  range  of  operating  frequency  over  which  the  half- power  beamwidth 
increases  no  more  than  25%,"  then  the  useful  bandwidth  can  be  stated  as 


Af  = 


f 

1.12  R 


(3-^D 
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If  the  half-power  sensitivity  bandwidth  is  used,  then 


Af  = 


O.56  R 


(3-42) 


One  should  note  that  the  present  "bandwidth"  considerations  are 
not  related  to  the  system  rise  time,  or  to  the  bandwidth  of  the  electro- 
acoustic sensor  channels  which  deliver  the  voltages  e*(qp,t)  in  Eq.  (3-l)- 
The  present  consideration  is  more  closely  related  to  correlator  mismatch 
error  between  the  assumed  and  operating  wavelength  of  the  signal  at  the 
array  face. 


The  effects  of  changes  in  the  operating  wavelength  are  more 
complicated  than  what  has  been  presented  here.  For  example,  if  the  shape 
of  the  major  lobe,  not  just  the  -3  dB  width  or  the  sensitivity  deterioration, 
is  considered,  then  a more  complicated  and  uncertain  prediction  of  behavior 
results.  This  effect  will  be  described  and  methods  of  reducing  the  effect 
will  be  considered  in  the  following  chapter. 

4.  Sensitivity — Array  Gain 

The  prediction  of  the  absolute  receiving  sensitivity  of  an 
array  (in  volts/plane  wave  p.bar)  is  very  difficult  because  of  its  dependence 
upon  several  factors,  such  as  the  absolute  sensitivity  of  the  elements,  or 
the  gain  of  the  preamplifiers.  Of  more  interest  here  is  the  relative  gain 
of  the  correlation  beamformer  when  compared  to,  say,  phase  shift  or  time 
delay  type  beamformers . 

In  general,  the  main  difference  between  the  correlation  method 
and  phase  shift  or  delay  line  methods  in  terms  of  sensitivity  is  that, 
where  the  two  latter  methods  move  each  element  electrically  so  that  all 
element  voltages  are  in  phase  and  sum  to  M (the  number  of  contributors) , 
the  correlation  method  adds  element  voltages  without  first  bringing  them 
into  phase  coherence . 
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If  the  collection  of  instantaneous  samples  of  element  voltages 
is  en>  then  the  signal  replica  factor  A^(n)  will  be  equal  to  e If  this 
effect  alone  is  considered,  then  the  output  voltage  would  be 

M 

E = £ A (n)  en  , (3-^3) 

n=l 


or,  since 


Ax(n) 


e 

n 


(3-44) 


E = 


(3-^5) 


Now,  if  the  array  is  large,  then  e^  is  a collection  of  samples  of  a 
cosine  wave  uniformly  distributed  throughout  the  cycle,  so  that  E might 
be  approximated  by 


or 


E 


2 

cos  cp  dcp 


E 


> 


(3-^) 


(3-^7) 


whereas,  the  phase  shift  or  time  delay  operation  results  in 


E = M . (3-48) 

The  general  result  then  is  that  the  sensitivity  of  an  array  with  a 
correlation  beamformer  is  at  least  6 dB  lower  than  coherent  processors, 
all  other  factors  being  equal.  Of  course,  given  a geometry  and  a beam 
pattern  specification,  the  other  shading  factors  will,  in  general,  be  the 
same. 

Another  item  of  interest  is  the  processing  gain  against  thermal 
noise.  Appendix  D contains  a detailed  comparison  between  a correlation 
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processor  and  a phase  shift  processor  which  were  designed  to  accomplish 
equivalent  beam  pattern  performance.  The  signal-to-thermal  noise  ratios 
were  within  0.5  dB  of  each  other. 

Methods  for  computing  sensitivity  and  S/N  ratio  are  outlined 
in  Appendix  D and  can  be  applied  to  individual  problems  as  the  need  arises . 


IV.  THE  DESIGN  AND  PERFORMANCE 
OF  MORE  REALISTIC  CYLINDRICAL  ARRAYS 


In  the  previous  chapter,  the  simplest  possible  cylindrical  array 
was  developed  in  order  to  convey  the  beamformer  concept,  to  show  the 
coefficient  selection  algorithms,  and  to  outline  the  performance  and 
scaling  rules.  This  array  was  called  the  reference  array. 

In  this  chapter,  we  will  consider  modifications  to  the  reference 
array  that  are  required  to  make  the  array  assumptions  more  realistic 
and  the  analysis  more  flexible. 

Figures  4-1  and  4-2  show  the  coordinate  system  used  in  this  and 
the  remaining  chapters  of  this  report.  Figure  3-1  will  not  be  repeated 
here,  but  is  applicable  as  well. 

A.  Finite  Size  Array  Sensors 

1.  Azimuthal  Effects 


One  of  the  first  modifications  to  the  reference  array  required 
to  make  the  analysis  more  realistic  is  the  inclusion  of  finite  size 
elements  or  staves.  Such  elements  exhibit  directional  response  indi- 
vidually, and,  since  they  are  assumed  to  be  mounted  on  the  cylindrical 
baffle  and  point  radially  outward,  the  effect  is  difficult  to  analyze. 
Intuition  indicated  that  the  element  beam  pattern  (azimuthal)  should  be 
accounted  for  in  the  direction  of  the  array  major  lobe,  since  such 
accounting  can  only  be  done  in  one  direction.  The  computer  model  was 
then  exercised  to  evaluate  the  array  pattern  performance. 
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For  the  farfield  case,  p and  g in  Fig.  4-1  are  very  large 
so  that  £n=cp-a^.  It  is  assumed  that  the  element  beam  pattern  is  symmetric 
about  its  pointing  direction  and  can  be  described  in  general  by 


B(?n)  = B(cp-an)  . (4-1) 

For  the  array  major  lobe  direction  (<p=0),  and  because  of 
the  element  pattern  symmetry, 


8(t„)  - B(an) 


(1.-2) 


Such  behavior  of  the  element  tends  to  shade  the  array  down 
toward  the  edges.  In  the  reference  case,  it  was  assumed  that  the 
elements  had  uniform  angular  response.  If  we  now  define  a new  shading 
factor  A(+(a),  called  the  "element  beam  pattern  corrector,"  and  further 
define  the  total  shading  coefficient  A(a)  as  (see  Eq.  J-2) 


A(a) 


A1(a)*A2(a)*A5(a) 


where 


Au(a)  = B(a)  = B(0 


(4-3) 


(4-4) 


then  we  will  have  raised  the  coefficient  value  at  an  angle  Ot in  such 
a k >y  as  to  account  for  that  element's  pattern  loss  in  the  direction  of 
the  array  major  lobe.  The  resulting  performance  should  be  nearly  the 
n*  at,  that  of  the  point  element  case  for  pattern  parameters  near  the 
major  lobe. 


Before  exercising  the  computer  model,  a specific  functional 
ff rm  for  the  B(a)  was  needed.  Prior  experience  with  small  width  elements 
(less  *han  0.5  \ wide)  on  rather  large  cylindrical  baffles  (R>10)  indicated 
that  the  pattern  could  be  adequately  described  by 
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B(d  = (X.  * f S g )’ 


where  x is  a constant  selected  to  make  B fit  a measured  element 
pattern. 


For  the  computer  runs  investigating  the  performance  of  finite 
size  elements,  several  values  of  x were  UBed.  Figure  4-3  shows  the 
beam  pattern  obtained  when  x=3,  and  should  be  compared  with  Fig.  3-4, 
the  point  element  counterpart.  The  sensitivity,  beamwidth,  and  peak  side 
lobe  levels  are  identical.  The  average  side  lobe  level  for  the  finite 
size  element  case  is  slightly  lower  because  of  the  pattern  smoothing 
evident  in  Fig.  4-3. 

In  summary  then,  one  can  design  arrays  using  the  point  element 
assumptions  and  scaling  rules.  During  construction  of  a real  array, 
element  patterns  can  be  measured  and  a functional  form  determined.  Then 
A^(a)  can  be  selected  and  incorporated  in  the  shading  coefficients  as 
shown  in  Eq.  (4-3)  without  affecting  the  point  element  performance 
prediction. 

2.  Elevation  Effects 

In  Chapter  III  we  presented  the  elevation  beam  pattern  of  the 
point  element  array  (Fig.  3-7)  and  the  scaling  rules  for  that  effect. 

V.hen  finite  size  elements  are  utilized  and  the  element  elevation 
pattern  is  B(*l)  (it  is  assumed  throughout  that  B(£,t|)=b(0*B(t)  )),  it 
was  found  that  the  composite  array  elevation  pattern  through  the  cp=0 
plane  was  very  closely  described  by  the  product  of  the  element  pattern 
B(*l)  and  the  point  element  array  pattern  (recall  the  line  array  product 
theorem) . 
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Figure  4-4  shows  a composite  elevation  pattern  with  element 
length  3.39  k • Such  an  element  alone  should  have  a sin  x/x  beamwidth  of 
15°  and  first  side  lobe  level  of  -13-5  dB.  The  comparable  point  element 
array  elevation  pattern  is  the  solid  curve  in  Fig.  3-7.  The  composite 
in  Fig.  4-4  has  a l4.8°  beamwidth  (not  much  affected)  and  a first  side 
lobe  of  -I6.5  (a  reduction  of  3 dB  due  to  the  array  response).  More 
will  be  said  about  this  effect  in  Chapter  VIII  where  experimental 
validation  of  the  theory  is  discussed. 

B.  The  Effei  ts  of  Element  Spacing 

In  order  to  determine  the  effects  of  various  element  spacings,  the 
computer  model  was  exercised  for  element  spacings  from  0.1  \ to  0.8 
holding  the  cylinder  radius  constant  at  R=12.3.  The  elements  were  points 
and  only  azimuthal  patterns  were  obtained.  Figure  4-5  shows  one  such 
pattern  with  a spacing  of  0.6  \. 

One  effect  that  can  be  seen  in  Fig.  4-5  is  the  growth  of  the  classic 
"grating  lobe"  at  about  120°.  As  the  spacing  grows  larger,  this  lobe 
becomes  larger  in  amplitude  and  moves  toward  the  major  lobe  at  0°.  The 
grating  lobe  behavior  is  identical  to  that  of  a phase  shift  or  time 
delay  beamformer. 

The  second  and  more  important  effect  seen  in  Fig.  4-5  is  the 
breakup  of  the  usual  uniform  side  lobe  region.  As  the  element  spacing 
increases  beyond  0.5  K,  the  peak  side  lobe  grows  rapidly.  A summary  of 
the  effect  is  shown  in  Fig.  4-6  for  an  R=12.3  array.  The  effects  shown 
here  indicate  that  when  side  lobe  level  is  important  (as  it  almost 
always  is)  one  must  space  elements  on  the  cylinder  at  less  than  O.55  \. 
and  preferably  at  0.5  \. 
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C.  Effects  of  A,  on  Beamwidth  and  Side  Lobe  Level 

Variations  in  azimuthal  beam  pattern  performance  of  the  correlation 
array  were  investigated  by  application  of  different  line  aperture  illumi- 
nation functions  (A,).  References  2 through  6 discuss  such  functions 
and  their  application  to  line  apertures.  Selected  amplitude  weighting 
of  the  element  outputs  produced  changes  in  beamwidth  and  peak  and  average 
side  lobe  level. 

The  study  involves  the  application  of  the  named  line  shading 
functions  to  the  equivalent  line  aperture  of  the  array  (the  W line).  The 
function  is  then  projected  back  to  the  cylinder,  where  the  A^  coefficient 
factor  is  evaluated.  The  computed  patterns  involve  the  total  composite 
coefficient  A{a)=k^{oi)*A^{a)*k^{OL)  as  defined  and  used  previously.  The 
A^  used  in  the  reference  array  was  a cosine  on  the  cylinder,  which 
results  in  a "house  top"  or  "triangle"  on  the  W line. 

The  effects  of  each  weighting  function  were  evaluated  by  a computer 
program  which  implements  a mathematical  model  of  the  cylindrical  array 
response  to  a point  source  in  the  farfield.  In  all  cases  the  array  had 
a diameter  of  100  \.  and  the  elements  were  separated  by  \/2,  thus 
resulting  in  628  elements  on  a l80°  aperture.  Following  is  a list  of 
the  line  array  weighting  functions  which  were  considered: 

1)  Uniform 

2)  Triangle  Functions 

3)  cos(np/2) 

4)  k + (l-k)  cos(nnp/2)  where  0 § k § 1 and  0.5  S n § 3 

Q 

5)  k + ( l-k)  cos^nnp ,f2)  (Hamming  Function)  where  0 5 k £ 1 and 
0.5  s n S 3 

6)  Taylor  Functions 

where  p is  the  normalized  distance  from  the  W line  center  to  the 
subtended  array  edge,  that  is,  0 £ p £ 1. 


*♦3 
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Figure  4-7  shows  the  results  of  applying  some  of  the  more  useful 
shading  functions.  The  bars  connect  the  peak  side  lobe  point  with  the 
corresponding  average  side  lobe  point.  Since  this  test  was  conducted 
on  an  R=100  array,  the  resulting  beamwidth  numbers  can  be  connected  to  K 
(the  aperture  constant)  by  multiplying  by  200.  The  point  to  be  emphasized 
in  Fig.  4-7  is  the  very  rapid  loss  in  the  aperture  constant  as  the  peak 
side  lobe  is  pushed  down. 

Equal  weighting  of  the  line  aperture  produces  the  narrowest  beamwidth 
(0.252°)  and  a very  low  average  side  lobe  level  (-26.2  dB),  but  with  a 
poor  peak  side  lobe  of  -14.56  dB.  Weighting  the  line  aperture  with  a 
simple  triangle  function  improves  the  response  to  a beamwidth  of  O.2960, 
the  average  side  lobe  level  of  -26.9  dB,  and  a peak  side  lobe  of  -20  dB 
(the  standard  A^=cos  a).  Other  often  used  line  weighting  functions,  such 
as  cosine  and  Hamming,  decrease  the  peak  side  lobes  to  -22.4  dB  and 
-23.4  dB,  but  with  a subsequent  growth  in  beamwidth  to  0.345°  and  O.380. 
respectively.  Of  the  functions  tested,  the  Taylor  functions  (Ref.  2) 
produce  the  best  compromise  between  beamwidth  and  side  lobe  level. 

D.  Effects  of  Changes  in  A^(a) 

The  reference  array  discussed  in  Chapter  III  assumed  the  signal 
to  be  positioned  on  the  array  so  that  a cosine  wave  existed  at  t=0 
(see  Fig.  3-3).  Then  A^(0f)  was  chosen  to  be  a replica  of  that  signal 
at  t=0.  Equation  (3-6)  suggests  that  other  times  could  be  chosen  for 
peak  response  by  using  the  cosine  phase  shift  ^ in  the  reference  signal. 

If  A^(ct)  is  then  matched  to  the  phase  shifted  signal,  it  becomes 

A^(a)  = cos  [2itR(i  - cos  a)  - . (4-6) 

Now  the  major  effect  of  changes  in  ^ is  the  shift  of  the  acoustic 
center  of  the  array;  however,  the  shift  is  not  in  direct  proportion  to 
changes  in  ♦.  The  explanation  can  be  found  in  the  shift  in  "center  of 
mass  of  the  coefficients." 
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FIGURE  4-7 

BEAMWIDTH  VERSUS  PEAK  AND  AVERAGE  SLL  ON  A CYLINDRICAL  ARRAY 


As  i‘  changes  from  zero,  the  original  large  coefficient  at  the  leading 
edge  of  the  array  (see  Fig.  3-M  is  reduced,  shifting  the  large  coeffi- 
cients farther  back  in  the  array.  The  beamwidth,  peak  side  lobe  level, 
the  overall  array  sensitivity,  and  the  acoustic  center  are  affected, 
although  for  the  first  three  this  is  a second  order  effect. 

These  effects  have  not  been  explained  analytically  and  graphs  of  the 
effects  are  not  shown  because  of  the  great  complexity  involved. 

Computer  designs,  however,  should  include  changes  in  ^ to  assess  the 
effect  on  the  beam  pattern  parameters  for  the  specific  case  at  hand. 

E.  Array  Focusing--The  Effects  of  Nearfield  Operation 

In  Chapter  III,  the  assumed  signal  originated  in  the  very  farfield, 
presenting  to  the  array  a plane  wave.  The  amplitude  of  the  signal  was 
uniform  over  the  width  and  depth  of  the  array,  and  the  phase  fronts  were 
plane.  In  this  section,  we  consider  the  generalization  of  the  coeffi- 
cient selection  algorithms  which  will  focus  the  array  at  a target  p 
(wavelength  units)  from  the  center  of  the  array.  The  nearfield  target 
will  be  assumed  to  lie  in  the  plane  of  the  array  (6=0°)  in  order  to 
simplify  the  description.  Figure  ^-8  shows  the  geometry  for  the  near- 
field target  including  the  curved  wavefront.  The  target  is  a point 
radiating  a monochromatic  wave  spreading  according  to  the  following 
expressions . 

P(g)  = — e^g  e^^  (linear  units)  , (4-7) 

o 

or 

P(G)  = e*^2nG  gJ^xcT  (wavelength  units)  . (4-8) 

0 
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The  coefficient  factor  A.,  (a  ) is  the  signal  replica  factor.  That  is 

1 n 

Ai(an)  = Real  Part  [Sig(cp=0,t=to)  ] (4-9) 

where  q*=0  is  the  beamforming  direction  and  t is  a time  selected  to 

cause  the  n=0  element  at  ao=0  to  see  a peak  positive  pressure.  That  is, 

t is  such  that 
o 

cos(2n[fto  - ( p-R) ])  = +1  . (4-10) 

The  amplitude  of  the  signal  is  also  normalized  to  unity  at  element 
n=0  by  the  factor  p-R. 

Now  the  range  to  any  element  can  be  found  by  the  law  of  cosines  as 


yp  _ p 

R + p 2Rp  cos  , (4-H) 

and  the  normalized  pressure  at  the  general  element  0^  is 

P(a  ) = cos [ 2n( ft  -G  )]  , ( 4-i 2) 

' n G on'  ' ' 

n 

or 

P(a  ) = cos[p-R-G  ] . (4-13) 

n ii  n 

n 

Equations  (4-12)  and  (4-13)  describe  the  normalized  signal  pressure 
at  each  element,  and  as  usual,  we  set 

W = P(°£n)  ' (4“lU) 

Now  A2(a)  is  the  sensor  density  corrector.  It  was  computed  for 
the  plane  wave  used  in  Chapter  III. A. 2.  The  same  algorithm  will  be 
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used  here,  but  with  different  geometry  as  shown  in  Fig.  4-8.  In  this 
case,  the  W line  is  curved  and  the  relationship  between  the  apparent 
element  position  on  the  W curve  (Y\  ) and  the  element  position  in  the 
array  (d^)  is 


Y.  = (p-R)Pi  , 


where  p can  be  found  by  the  law  of  sines,  resulting  in 


B.  = arcs in 
i 


[£--  sin  0./j  , 


or 


Yi  = (p-R)  arcsin  sin  <xj 


(4-15) 


(4-16) 


(4-17) 


Now,  following  the  outline  in  Chapter  III. A. 2.,  if  the  element 
density  on  the  cylinder  is  q/d (X,  then  the  apparent  density  on  the  W 
curve  is  q/dY.  The  derivative  of  Eq.  (4-17)  with  respect  to  a , 
remembering  that  G is  a function  of  0!  also,  will  allow  the  evaluation  of 
q/dY  in  terms  of  q/d a.  In  addition  to  this  density  term,  the  envelope 
of  both  the  signal  and  A^(a)  should  be  included.  The  normalized 
signal  strength  at  the  nth  element  is 


£_ 1 


R 


G 


n 


p + R 


(4-18) 


and  A^  was  chosen  to  have  the  same  envelope  (Eq.  4-l4).  So,  although 
the  apparent  element  density  near  the  edge  of  the  W curve  increases, 
the  signal  and  A^  envelopes  decrease.  Thus,  A0(a)  should  be  chosen  to 
make 

Ap(a)  f — P ^ = constant  . (4-19) 

lGn  ' P + RJ 
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Now,  Aj(a)  is  the  aperture  window  function  and  should  be  applied 
to  the  W curve.  Its  value  at  each  element  can  be  evaluated  using 
Eq.  (4-17)  to  relate  a to  Y. 

Also,  Aj^a)  (the  element  beam  pattern  corrector)  should  be  modified 
slightly  to  include  the  effect  of  parallax.  As  can  be  seen  in  Fig.  4-8, 
the  correction  to  the  element  pointing  angle  is  simply  6^,  so  that 
Ait(a+p)  is  the  element  beam  pattern  corrector  where  (3  is  given  in 
Eq.  (4-16). 

Experiments  and  computer  modeling  indicated  that,  if  the 
coefficient  selection  algorithms  described  above  are  utilized,  then 
the  performance  of  a focused  beam  with  a target  at  the  focused  range 
is  essentially  the  same  as  the  performance  of  the  unfocused  reference 
beam  of  Chapter  III  with  a target  in  the  farfield. 

F.  Frequency  Response 

In  Chapter  III.B.3-,  the  bandwidth  of  the  reference  array  was 
discussed.  The  array  sensitivity,  half-power  beamwidth,  and  the  peak 
and  average  side  lobe  changes  due  to  changes  in  operating  frequency 
were  shown.  A more  stringent  requirement  on  the  operating  frequency 
may  be  the  change  in  shape,  or  in  form  factor,  of  the  major  lobe. 

In  order  to  study  this  effect,  a large  array  (R=136)  was  chosen. 
Figure  4-9  is  a plot  of  the  beamwidths  measured  at  different  levels  on 
the  major  lobe  as  a function  of  operating  frequency.  The  array  was 
designed  to  operate  at  224  kHz.  The  sector  of  elements  covers 
23^=120°,  so  the  array  depth  (R( 1-cos  0^)= 68)  is  68  wavelengths. 

Eq.  (3-4l)  would  indicate  a 2.94  kHz  operating  band  in  which  the  -3  dB 
beamwidth  would  not  grow  beyond  25$.  This  estimate  is  confirmed  in 
Fig.  4-9.  However,  over  this  band,  the  -15  dB  beamwidth  is  seen  to 
grow  by  about  45$.  in  some  cases,  it  may  be  necessary  to  restrict  the 
operating  bandwidth  further  to  control  the  major  lobe  shape.  If  the 


operating  band  is  restricted  to  one-half  that  predicted  by  Eq.  (3-4l) 
(~1.5  kHz),  then  even  the  -15  dB  beamwidth  grows  less  than  20 4o. 

The  three  arrows  shown  in  Fig.  4-9  indicate  the  three  operating 
frequencies  used  for  the  beam  patterns  shown  in  Figs.  4-10,  4-11,  and 
4-12.  It  can  be  seen  that  the  peak  side  lobe  varies  somewhat  over  the 
band,  but  the  average  side  lobe  is  constant.  The  -3  dB  beamwidths  are 
listed  on  each  figure. 

G.  The  Effects  of  Errors 


This  section  concerns  the  effects  on  the  beam  pattern  parameters 
caused  by  either  design  or  construction  errors  in  either  the  element 
phase  and  amplitude  response,  element  placement  errors,  or  the  shading 
coefficients . 

Equation  (3-l)  represents  the  perfect  beamforming  operation. 

It  is  repeated  here  as 

* N * 

Eo  (<P»t)  = £ An  en  (*>*)  > (4-20) 

n=-N 

where  * denotes  a complex  number.  If 

e (qp, t ) = e exp[i0  ] , (4-21) 

n ’ n n ’ 

and  all  possible  amplitude  errors  are  grouped  in  6^,  ana  all  phase 
errors  are  grouped  in  Ar,  then  the  imperfect  beamformer  operation  can 
be  written  as 


L 

n 


e 6 
n n 


exp[i(B  -A  ) ] 
n n 


(4-22) 


52 


0 


-10 


CD 

-a 


>■ 

h- 

> 


«/> 
Z 
LU 
i n 


>■ 

< 

on 

on 

< 


-20 


-30 


-40 


u 

0 2 

4 6 

ANGLE  - deg 

8 

10 

RADIUS  = 136.55  A 

BEAMWIDTH  * 

0.2962* 

RADIUS  = 91.44  cm 

PEAK  SLL  * 

-23.0  dB 

FREQUENCY  (DESIGN)  = 
FREQUENCY  (OPERATE) 

224.0  kHz 
= 223.5  kHz 

AVERAGE  SLL 

: - 25.3  dB 

FIGURE  4-10 

BEAM  PATTERN  AT  223.5  kHz 
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BEAM  PATTERN  AT  224  kHz 
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BEAM  PATTERN  AT  224.5  kHz 
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The  first  test  of  the  effects  of  error  was  performed  with  uniform 

distributions  for  the  errors.  In  that  test,  6 was  a random  variable 

n 

with  uniform  distribution  whose  mean  was  unity  and  whose  extreme  values 
were  1±X  where  100  X was  the  percentage  amplitude  error  under  considera- 
tion. The  random  variable  A was  also  drawn  from  a uniform  distribution 

n 

with  zero  mean  and  extreme  values  ±Y  stated  as  phase  error  in  electrical 
degrees . 

Figure  4-13  shows  an  azimuthal  beam  pattern  with  6^  drawn  from  a 
uniform  distribution  from  0.6  to  1.4,  called  the  ±40$  error,  and  A sO. 

The  figure  should  be  compared  to  Fig.  3-5>  which  is  the  zero  error  case. 
The  beamwidth  is  slightly  larger  in  Fig.  4-13,  the  sensitivity  is  slightly 
greater,  and  the  side  Lobe  region  is  more  scattered  than  in  Fig.  3-5- 

Figure  4-14  is  an  example  summary  plot  showing  the  effect  of 
j varying  the  magnitude  of  the  amplitude  error.  Only  five  beam  patterns 

were  computed  with  each  error  value.  Other  summary  plots  of  the  type 
shown  in  Fig.  4-14  with  various  amplitude  and  phase  errors  were 
generated,  with  the  result  that  a ±40°  phase  error  had  generally  the 
same  effect  on  SLL  and  SLL^  as  a ±60$  amplitude  erz’or. 

The  array  sensitivity,  on  the  average,  will  depend  upon  the  average 
sensitivity  of  the  elements  and  the  number  of  contributors.  Using  the 
same  arguments,  the  average  side  lobe  level  depends  upon  the  incoherent 
addition  of  the  element’s  contributions.  The  variance  of  the  side  lobe 
level  grows  with  element  amplitude  error  but  decreases  with  the  number 
of  contributors.  The  average  side  lobe  level  decreases  with  an  increase 
in  the  number  of  contributors.  To  the  first  order,  then,  the  variance 

I of  the  average  side  lobe  level  depends  only  upon  the  element  error,  not 

upon  the  size,  number  of  elements,  or  the  absolute  side  lobe  level. 

This  statement  is  demonstrated  in  Fig.  4-15,  which  shows  two  patterns 
with  different  numbers  of  elements  (and  hence  different  average  side 
lobe  levels)  but  with  the  same  element  error.  One  can  see  the  side  lobe 
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AZIMUTHAL  BEAM  PATTERN 
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variance  is  nearly  constant  in  the  decibel  scale.  The  conclusion  that 
can  be  drawn  from  the  above  argument  is  that  larger  arrays  with  lower 
side  lobe  levels  do  not  require  the  use  of  a more  uniform  set  of 
elements . 

A second  study  of  the  effects  of  error  is  summarized  in  Figs.  ^-l6, 
4-17,  and  4-l8.  In  this  case,  a larger  array  was  used  (R=l20,  20^=120°, 
containing  5^3  elements).  The  errors  8^  and  were  taken  from  a normal 
distribution  whose  standard  deviation  is  stated  as  the  error.  In  each 
case  tested,  25  beam  patterns  were  computed  so  that  averages  and  standard 
deviations  of  the  effects  were  plotted  with  moderate  confidence. 

Figure  4-16  shows  a plot  of  several  parameters  versus  amplitude 
error.  Figure  ^-17  shows  the  effects  of  phase  error,  and  Fig.  ^-18 
combines  the  effects  by  assuming  a 25$  amplitude  error  (a=0.25)  and 
plotting  the  additional  effects  of  increasing  phase  error.  It  can  be 
seen  that,  for  25$  amplitude  and  10°  phase  error  (a=10°),  the  peak 
side  lobe  of  80$  of  the  computed  patterns  fell  within  1.5  dB  of  the 
pattern  with  no  error. 
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V.  THE  ARRAY  FILL  TIME  PROBLEM  AND  PROPOSED  SOLUTIONS 


A.  The  Fill  Time  Problem 


If  an  array  of  sensors  used  for  beamforming  subtends  some  nonzero 
dimension  in  the  direction  of  signal  propagation,  and  does  not  utilize 
time  delay  processing,  then  it  will  exhibit  some  array  rise  time.  The 
effect  is  caused  by  the  time  required  for  the  signal  to  fully  in6onify 
the  array.  During  the  rise  time,  the  sensitivity  is  lower  than  usual 
(because  not  all  elements  are  insonified),  and  the  beamwidth  is  larger 
than  usual  (because  only  the  center  section  of  elements  making  up  a 
smaller  aperture  are  insonified).  This  effect  is  separate  and  distinct 
from  the  effect  of  signal  rise  time  caused  by  bandwidth  restrictions  in 
the  receiver  channels. 

The  solid  curve  in  Fig.  5-1  shows  the  speed  with  which  the  array 
output  level  rises  as  a rectangular  pulse  insouifies  the  reference  array. 
The  rapid  rise  and  fall  of  the  output  level  is  due  mainly  to  the  weight 
of  array  elements  near  the  front  center  which  are  insonified  first  and 
which  are  first  to  be  uninsonified  as  the  pulse  leaves  the  array.  The 
dashed  curve  in  the  figure  shows  the  growth  in  the  aperture  W as  the  pulse 
comes  on  the  array.  The  resultant  relative  beamwidth  is  related  to  the 
reciprocal  of  the  aperture. 

Another  point  of  interest  in  Fig.  5-1  is  that  the  apparent  pulse 
length  is  very  nearly  equal  to  the  real  pulse  length,  unelongated  by  the 
array.  The  conclusion  is  that  the  range  resolution  achievable  with  the 
simple  cylindrical  array  is,  as  with  other  arrays,  related  to  one-half 
the  pulse  length. 
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If  the  pulse  just  covers  the  array  radius  (the  minimum  usable  length), 
the  pulse  length  t is 


T = £ _ R 
T " c " f 


(5-D 


and  then  the  range  resolution 


A CT 

Ax  = — 


(5-2) 


is 


. c R r 
A X = 2 f = 2 


(5-3) 


Now,  the  crossrange  resolution  Ay  at  a target  range  p yards  from  a 
sonar  with  beamwidth  BW  degrees  is  given  by 


Ay  = pBW  yards 


(5-4) 


But,  for  the  correlation  beamformer, 


BW  = ~ degrees  , 


(5-5) 


where  K is  the  aperture  constant.  This  results  in  the  expression 


pIQt 


^ - ^Tr  * 5 feet  ' 


(5-6) 


for  the  crossrange  resolution.  We  can  now  compute  the  area  of  a 
resolution  cell  at  range  p yd  as 


Area  = AxAy  = ttc 


(5-7) 


Area  = ^ x( constant) 


(5-8) 
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The  point  is  that  the  area  of  a resolution  cell  is  independent  of  the 
array  size.  That  is,  as  one  constructs  larger  arrays  to  achieve  better 
crossrange  resolution,  the  pulse  length  must  grow  to  cover  the  larger 
array,  with  the  result  that  the  resolution  area  at  some  range  p remains 
constant . 

In  summary,  if  the  sonar  designer  wishes  a system  with  pulse  length 
PL  (seconds)  and  beamwidth  BW  (degrees)  operating  at  frequency  f (hertz), 
and  if  the  product 

f * PL  * BW  S 30  , (5-9) 

he  can  use  the  correlation  beamformer  with  a simple  rectangular  CW  pulse 
and  2a^=l80°.  The  array  diameter  will  be 

2Kc 

2r  = feet  . (5-10) 

o 

If,  on  the  other  hand, 

f * PL  * BW  < 30  , (5-11) 

then  the  required  pulse  length  will  be  too  short  for  the  array.  This 
represents  the  rise  time  problem  for  the  correlation  array,  and  two 
solutions  are  presently  recognized  as  possible.  They  are  (l)  reduction 
of  the  depth  of  the  array  by  reducing  20^  or,  (2)  utilizing  a coded 
pulse  transmission  and  performing  a replica  correlation  on  the  receiving 
array.  The  two  methods  will  be  treated  separately  in  the  following  two 
sections. 


B.  The  First  Solution — Reduced  Sector  of  Elements  in  the  Aperture 


One  suggested  solution  to  the  rise  time  or  range  resolution  problem 
i s to  reduce  the  sector  of  active  elements  on  the  cylinder.  Normally, 
20^=180°,  which  means  that  a pulse  must  be  long  enough  to  insonify  the 


a 
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radius  of  the  array  for  on-axis  targets.  If  the  sector  angle  20^  were 
reduced,  the  pulse  length  required  to  insonify  the  array  would  be 
reduced.  However,  a disadvantage  of  this  procedure  is  the  attendant 
reduction  in  the  effective  beamforming  aperture  W if  R is  not  increased 
with  a reduction  in  L. 

The  depth  of  aperture  is  given  by 

L = R(1  - cos  otjj)  . (5-12) 

The  effective  beamforming  aperture  W is  given  by 

W = 2R  sin  aN  , (5-13) 

where  a„  is  one-half  the  total  angular  sector  of  elements  on  the  array. 

N 

Equations  (5-12)  and  (5-13)  can  be  rewritten  to  show  the  ratios 

I = 1 - cos  aN  (5-iM 

and 

7$  = sin  aN  , (5-15) 


and  to  illustrate  how  these  ratios  depend  upon  0^. 

Now,  if  one  wishes  to  decrease  L in  order  to  shorten  the  apparent 
pulse  length  while  preserving  the  beamwidth,  then  one  must  solve  both 
Eqs.  (5-12)  and  (5-13)  simultaneously.  Solving  Eq.  (5-13)  for  R yields 


R = 


W 

2 sin  0^ 


Substituting  Eq.  (5-16)  Into  Eq.  (5-12)  yields 

T wA  - COB  °^\ 

= 2\  / 


(5-16) 


(5-17) 
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which  illustrates  the  reduction  in  L achievable  by  reducing  a^. 

Figure  A-4  in  Appendix  A is  a plot  of  the  angle  function  expressed  in 
Eq.  (5-17).  By  entering  that  plot  with  the  required  fractional  reduc- 
tions in  L below  the  original  R value,  the  corresponding  angle  can 
be  obtained.  This  angle  can  then  be  substituted  into  Eq.  (5-16)  to  find 
the  new  value  of  R. 


Small  reductions  in  L below  the  2a^=l80°  value  can  be  accomplished 
with  little  increase  in  R;  however,  further  reductions  in  L lead  finally 
to  alarming  increases  in  R. 

This  method  of  resolving  the  array  rise  time  problem  is  rather 
limited  and  is  costly  in  other  areas  as  well  as  those  discussed  above. 

For  example,  there  is  an  increase  in  the  number  cf  elements/ staves  and 
associated  preamplifier  channels  required  when  several  preformed  beams 
are  considered.  Appendix  A contains  a discussion  of  these  areas  and  other 
associated  factors  when  the  sector  of  active  elements  is  reduced,  for 
whatever  reason. 


C.  The  Second  Solution — A Coded  Pulse  Transmission  and  Replica 

Correlation  at  the  Receiving  Array 

Because  the  array  has  depth  in  the  direction  of  signal  propagation, 
it  can  be  used  as  a correlator.  Correlation  is  normally  imnlemented  with 
some  sort  of  delay  element  (delay  line,  shift  register,  etc.)  and  a 
summer.  We  recognize  that  the  delay  element  in  our  case  is  the  array 
itself,  with  the  signal  propagating  past  the  array  at  the  speed  of  sound. 
Since  the  array  has  depth  L,  it  can  simultaneously  sample  as  many  as  L 
cycles  of  the  signal.  In  the  reference  array  case,  20^=180°,  and  L=R. 

In  Fig.  5-2,  we  imagine  the  array  as  being  collapsed  on  the  x-axis, 
occupying  the  region  from  the  origin  (the  array  center)  out  to  R (in 
wavelengths).  The  leading  edge  of  the  pulsed  signal  begins  at  g distance 
from  the  origin  at  time  t.  The  signal  description  is  S(r)  in  a moving 
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coordinate  system,  or  S(t+r)  in  the  array  centered  system.  Imagine  t 
time  running  down,  so  that  the  signal  approaches  the  array  from  the  right. 
When  t=T,  the  leading  edge  of  the  signal  just  touches  the  array.  During 
the  time  -T<t<T,  the  signal  partially  overlaps  the  array. 


In  Chapter  III,  we  argued  that  the  beamformer  performs  a multiply 
operation  on  the  sampled  signal  and  that  these  scaled  outputs  are  summed. 
In  the  limit  of  a large  number  of  elements  or  samples  then,  the  moving 
signal  passing  the  array  can  be  shown  to  cause  an  output  from  the  beam- 
former  of 


E(t) 


A(t)  S(t+f)  dr 


(5-18) 


Now,  if  we  wish  the  array  output  to  be,  for  example,  the 
self-convolution  of  the  signal,  then  we  could  choose  A(t)=S(t),  (the 
signal  replica)  and 


E(t) 


S(t)  S(t+x)  dT 


(5-19) 


In  this  case,  pulse  compression  could  be  accomplished  by  choosing  the 
signal  S(t)  to  be  any  code  with  good  autocorrelation.  Of  course, 

Eq.  (5-l8)  is  only  true  for  signals  arriving  on  the  MRA,  since  that  is 
the  only  arrival  angle  where  the  coefficients  A(t)  are  chosen  to  match 
the  signal. 


Now,  several  complications  arise  in  practice  because  the  angular 
response  of  the  array  must  also  be  considered.  It  is  conjectured  that 
if  the  coded  pulse  is  fully  incident  on  the  array,  that  is,  t=0,  and  the 
pulse  is  phase  coded  digitally  so  tb  .0  the  carrier  wavelength  is  constant 
throughout  the  pulse,  the  resultant  fully  insonified  azimuthal  beam  pattern 
will  not  be  too  differed':  from  the  long  CW  signal  case.  This  conjecture 
allows  us  to  separate  the  time  and  angle  response  of  the  array  during  the 
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design  phase,  and  then  to  test  the  conjecture  by  computer  trial.  An 
example  will  be  presented  later. 


In  view  of  the  conjecture,  one  can  first  design  the  azimuthal 
pattern  with  a long  CW  pulse  and  determine  the  shading  coefficients  as 
described  in  Chapters  III  and  IV  [they  are  A(an))-  These  coefficients 
expressed  in  the  time  domain  are  A(t)  and  are  seen  to  be  quite  different 
from  S(t).  This,  then,  is  the  major  problem  with  choosing  a suitably 
coded  pulse;  that  is,  Eq.  (5-18)  becomes 


where 


E(t) 


A'(t)  S(t)  S(t+x)  dx 


(5-20) 


A' (r)  S(t)  = A(x) 


(5-21) 


Equation  (5-20)  is  no  longer  a simple  self- convolution,  like  Eq.  (5-19) > 
but  a complicated  crosscorrelation  with  a highly  asymmetric  window 
function,  A'(x).  The  object  then  is  to  find  a phase  coded  pulse  S(x) 
such  that  the  E(t)  has  the  appropriate  pulsewidth  with  low  time  side  lobes. 


At  this  point  in  the  study,  no  organized  method  of  determining  S(x) 
has  been  found.  A trial  and  error  approach,  modified  somewhat  by 
intuition,  has  led  to  a code  class  that  shows  promise,  but  is  surely  not 
optimum.  The  following  paragraphs  outline  the  code  selection  procedure, 
and  the  results  of  computer  trials  are  shown  in  the  form  of  contour  maps 
of  E(cp,t),  showing  the  composite  time  and  angle  response  of  the  array 
with  various  signal  types.  The  first  section  will  consider  an  uncoded, 
short,  CW  pulse  with  rectangular  envelope.  The  major  difficulty  with 
rectangular  envelopes  will  be  described.  The  second  section  will  consider 
a phase  coded,  triangular  envelope,  CW  pulse,  which  exhibits  the  desired 
pulse  comoression. 
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1. 


The  Rectangular  Envelope  CW  Pulse 


Consider  a signal  S(t)  with  rectangular  envelope  of  length 
(T/2)+e  sec,  which  is  to  cover  the  front  120°  of  elements  on  a cylindrical 
array  of  radius  R.  The  signal  envelope 

|s(t) | = 1 when  J § r S T + e , (5-22) 

where  e will  be  explained  later,  and 

|S(t) | = 0 otherwise.  (5-25) 

When  the  signal  is  full  on  the  array,  t=0  and,  aligned  to  <po=0,  the  signal 
replica  factor  A^(a)  can  be  found  as  shown  in  Chapter  III.  As  can  be  seen 
from  Pig.  5-2,  when  t=0. 


which  leads  to 


and 


r cos  a = ct 


r = T cos  a 


a = arccos 


(5-2*0 


(5-25) 


(5-26) 


These  relations  allow  the  change  of  coordinates  required  to  change 

Ax(a)  = |s(a)|  cos[2rtR(l  - cos  a)]  (5-27) 

to 

A1(t)  = |S(t) | cos[2rtR(l  - i)]  . (5-28) 


Following  the  development  shown  in  Chapter  III. A. 2.,  the  sensor 
density  corrector  can  be  found.  Equation  (5-10)  shows  the  effect  of  the 
element  density  alone.  Actually,  we  should  include  the  envelope  of  the 
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signal  and  in  the  sensor  density  corrector.  Thus 

|A1(Y) | |S(Y)|  ± 

represents  the  effective  A^*  signal  aperture  function.  Of  course,  in 
this  case 

|a1(t)|  = |s(x) | = 1 

so  that,  just  as  in  Chapter  III, 

A2(a)  = cos  a 


since  it  causes 


|A1(Y) I |S(Y) I ^ Ag(Y)  = constant. 

A^(Y)  is  the  aperture  window  function  and,  in  the  reference  array, 
was  chosen  to  be  a "house  top,"  \J  1-(Y/R)2,  which  leads  to  A^(a)=cos  a. 

In  our  present  case,  we  let  A^(Y)  be  a Taylor  shading  function  with  a 
design  side  lobe  level  of  27  dB  and  N=l4.  This  function  is  reflected 
back  to  the  cylinder  and  sampled  to  obtain  A^(a).  The  Taylor  function 
and  its  reflection  to  the  cylinder  are  too  complicated  to  present  here; 
however,  the  method  should  be  clear. 

Now,  the  considered  array  is  R=320  and  2aN=120o,  with  503  point 
elements  spaced  \/2  apart,  the  signal  is  a cosinusoid  with  rectangular 
envelope  ac  described  previously,  and  the  coefficients  sire  selected  as 
described  above.  The  resulting  azimuthal  beam  pattern  is  shown  in 
Fig.  5-3-  It  will  be  remembered  that  the  pulse  is  e sec  longer  than  the 
depth  of  active  elements.  In  this  case,  the  pulse  was  5 cycles  longer 
than  the  array  depth  (60  \) , giving  a total  pulse  time  duration  of  325  psec. 
The  stated  time  in  the  figure  is  time  from  first  signal  insonifi cation  on 
the  array  and  is  one  pulse  length  different  from  the  t time  in  Fig.  5-2. 
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Therefore,  the  time  indicated  on  the  plots  V is  related  to  t time  by 


r = | - t . (5.29) 

Also,  as  t time  runs  down,  T time  increases.  Of  course,  full  insonifica- 
tion  occurs  during  the  time  duration  e from 

0 * t § - e , (5-30) 

or 

|srs|+e  = 325  nsec  (5-3l) 


for  this  case. 

Figure  5-3  thus  represents  a fully  insonified  beam  pattern. 

The  3 dB  beamwidth  is  less  than  0.3°  (O.29880),  and  the  minimum  side  lobe 
suppression  is  21.27  dB. 

Figure  5-^  shows  the  time  response  of  the  uncoded  pulse  case 
for  the  signal  arriving  on  the  major  lobe,  or  from  the  cpQ=0  direction. 

The  plot  starts  atV=25  nsec,  where  the  signal  is  already  5 cycles  into 
the  array,  and  at  that  time  the  array  response  is  already  within  10  dB 
of  the  peak  response. 

Questions  such  as  "What  is  the  beamwidth  of  the  azimuthal 
pattern  when  the  pulse  is  partially  on  the  array?"  or  "What  is  the  time 
duration  of  some  particular  side  lobe?"  can  be  answered  by  combining  beam 
patterns  measured  at  many  different  T times  into  a contour  map  of  the 
array  response  as  a function  of  V time  and  cp  angle.  Figure  5-5  is  such 
a contour  map.  It  is  instructive  to  compare  Figs.  5-3  and  5-4  with 
Fig.  5-5,  inasmuch  as  they  are  simply  cuts  from  the  contour  map  at 
T=325  usee  and  at  <p=0.  All  response  parameters  of  interest  can  be  read 
from  a plot  of  the  type  shown  in  Fig.  5-5* 
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So  far  in  this  chapter  no  attempt  has  yet  been  made  to  use  pulse 
jompression.  The  previous  paragraphs  serve  only  to  introduce  notation  and 
plotting  format,  and  to  show  the  resulting  response  if  no  pulse  compression 
is  utilized.  If  a shorter  effective  pulse  length  is  required  at  the  same 
beamwidth  and,  therefore,  array  size,  then  pulse  compression  will  be 
necessary.  In  spite  of  the  possible  existence  of  many  code  classes  which 
will  achieve  useful  pulse  compression  on  a cylindrical  array,  only  one 
has  so  far  been  found  . It  will  be  described  in  the  following  paragraphs . 

2.  The  Triangular  Envelope — Phase  Coded  Pulse 

The  object  of  pulse  compression  is  twofold.  First,  the 
apparent  pulse  length  out  of  the  beamformer  is  required  to  be  shorter 
than  the  real  pulse  length  in  the  water.  Second,  it  is  required  that  at 
all  other  times  (other  than  the  peak)  the  response  be  less  than  some 
level  (the  time  side  lobe  level). 

When  phase  coded  pulses  were  first  investigated,  it  was  found 
that  the  first  time  side  lobe  could  not  be  adequately  controlled  by  phase 
changes.  Apparently  the  reason  was  that  a rectangular  pulse  (large 
amplitude  at  the  beginning)  and  the  array  with  large  weights  concentrated 
near  the  leading  edge  combined  to  cause  a large  output  as  the  pulse 
entered  the  array.  Of  course,  we  would  like  the  array  response  to  remain 
low  until  the  pulse  covers  the  array.  This  early  response  peak  or  time 
side  lobe  can  be  controlled  by  utilizing  a triangular  envelope  signal. 

The  slowly  increasing  signal  strength  gives  sufficient  time  for  phase 
changes  to  be  effective  in  controlling  the  resulting  time  side  lobes. 

Given  a triangular  envelope  for  the  signal,  the  next  step  was 
to  determine  the  locations  of  l8o°  carrier  phase  shifts  that  would  keep 
time  side  lobes  within  bounds  and  allow  the  peak  response  to  remain  up 
for  a useful  time.  This  selection  was  accomplished  empirically,  essen- 
tially a trial  and  error  approach.  If  only  two  phase  changes  are  allowed, 
and  the  code  is  represented  by  1 3 equal  time  increments  or  bits,  then  the 
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best  code  found  is 


++++ — +++++++  , (5-32) 

which  shows  that  the  phase  changes  occur  after  the  fourth  and  sixth  bits. 

If  we  apply  such  a coded  pulse  to  the  previous  case  of  R=120, 
2aN=120° , and  503  elements  spaced  \/2  apart  at  fQ=200  kHz,  we  find  that 
the  active  depth  of  array  is  R cos  = 60.  If  we  let  e be  the  time 
required  for  5 extra  cycles,  the  total  pulse  length  can  be  broken  into 
13  equal  time  segments  or  bits  of  5 cycles  each.  The  duration  of  each 
bit  (5  cycles)  sets  the  transmitter  and  receiver  bandwidth  required  to 
reproduce  the  signal. 

Figure  5 -6  shows  a time  plot  (upper)  of  the  unfiltered  phase 
coded  pulse  of  325  nsec  length.  Just  below  is  shown  the  time  plot  of  a 
filtered  version  (simulating  the  composite  transmitter  and  receiver 
bandwidth  of  32  kHz)  of  the  same  coded  pulse.  It  is  slightly  delayed  and 
elongated  by  the  filter.  Let  S(t')  represent  the  unfiltered  signal  and 
S'(t')  represent  the  filtered  waveform.  The  shading  coefficients  will 
now  be  computed  and  the  results  shown. 

The  signal  replica  A^(x)  can  be  found  by  setting 

A1(t)=S’(t')  (5-33) 

where  the  relationship  between  t and  x'  is  most  important  and  is  responsible 
for  the  signal  placement  on  the  array  at  the  time  of  maximum  response. 

Let  us  choose  ^ 

r'  = x - | + € , (5-3*0 

making  t=0  for  the  time  of  maximum  response,  as  shown  in  Fig.  5-7*  In 
the  case  under  consideration 
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FIGURE  5-6 

CODED  SIGNAL  PULSE  WAVEFORM 
BEFORE  AND  AFTER  FILTERING 


FIGURE  5-7 

REPLICA-CORRELATION  GEOMETRY 
FOR  CODED  WAVEFORM 
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2 = 500  Msec  , (5-55) 

and  e is  the  time  required  for  5 cycles.  Therefore, 

e = 25  qsec  (5-36) 

and 

r'  = t - 275  dsec  . (5-37) 

The  upper  scale  in  the  lower  plot  of  Fig.  5-6  illustrates  the 
relationship  between  x and  x1  at  t=0,  the  time  of  maximum  response.  Also, 
the  arrows  at  x=500  and  x=600  show  the  portion  of  the  signal  incident  on 
the  array  at  t=0. 


The  array  is  cylindrical  with  equally  spaced  elements,  so  that 
the  element  density  along  the  W line  is,  as  before, 

3y  = da  R cos  a * (5-38) 

But  this  time,  the  signal  envelope  and  A^  envelope  are  not  constant  in 
time.  We  should  use  the  detailed  envelope  of  the  filtered  signal  S'(x'); 
however,  because  of  its  complexity,  we  will  say  that  the  signal  has  an 
average  envelope  like  S(x'),  which  is 

I ^(t  1 ) I — ™ — , (5-39) 

— + e 


when 


0 S x ' 5 i + € 


By  using  Eqs.  (5-33)  and  (5-34),  Eq.  (5-39)  becomes 


2x  - T + 2e 


2e 


(5-40) 


(5-41) 
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where 


T 

i S t S T 


(5-te) 


J 


If  we  divide  by  T in  Eq.  (5-^l)  and  use  the  value  of  T and  e except  in  the 
first  term  where  t/T=cos  a is  substituted,  we  obtain 


|A1(a)|  = 


2 cos  a - 

n 

12 


(5-^3) 


or 

lAl(a)  | = jj  (cos  a - , (5-W) 

where 

0 s a s 6o°  . (5-^5) 

Of  course,  the  signal  description  on  t in  the  same  range  is  the  same. 


Now  we  have  enough  information  to  determine  so  that 

|A1(Y)  I |S(Y)|  Ag(Y)  = constant.  (5-^) 

From  Eq.  (5-^),  ignoring  the  constant  factor,  and  Eq.  (5-38),  again 
ignoring  the  constant  factor,  we  have 


Ag(a)  = 


Q cos  a 


fos  « - 


(5-^7) 


where  Eq.  (5-^*6)  has  been  couched  in  terms  of  a and  where  Q is  the  constant 
in  Eq.  (5-^6),  but  which  we  will  use  to  normalize  Ag  at  a«0.  Thus, 


(5-W) 
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Now  A^(a)  will  be  chosen,  as  in  the  example  without  pulse 
coding,  to  be  a Taylor  function  with  design  side  lobe  level  of  27  dB 
and  N=l4  on  the  W line,  and  then  reflected  to  the  cylinder. 

The  total  coefficients  can  now  be  computed  by 

A (a)  = A1(a)  * A^a)  * A?(a)  . (5-49) 

Figure  5-8  is  the  resulting  azimuthal  beam  pattern  with  pulse 
coding  at  t=0  or  V=325  usee.  The  beamwidth  is  0.3235°  with  the  largest 
side  lobe  suppressed  to  a level  21.87  dB  below  the  peak.  It  appears  that 
the  beamwidth  grew  about  8$  over  the  uncoded  case.  The  side  lobe  level, 
however,  is  slightly  lower. 

Figure  5-9  shows  the  time  response  for  the  coded  pulse  case 

at  <p=0.  The  two  early  time  side  lobeB  are  suppressed  more  than  15  dB, 

% 

and  the  -6  dB  pulse  length  measures  120  psec.  The  large  back  time  side 
lobe  is  caused  by  the  choice  of  t time  of  maximum  response,  of  which 
more  will  be  said  later. 

Finally,  Fig.  5-10  is  the  contour  map  of  the  beamformer  output 
plotted  against  both  time  and  angle.  It  should  be  compared  to  Fig.  5~5» 
which  shows  the  uncoded  case.  It  should  be  noted  that  Fig.  5-10  does 
not  represent  the  best  pattern  achievable;  however,  it  is  a pattern  of 
some  use  and  has  the  distinct  advantage  of  having  an  explainable  lineage. 

A summary  of  the  choices  made  in  order  to  arrive  at  Fig.  5-10 
may  convey  the  number  of  parameters  that  need  to  be  investigated  in  order 
to  arrive  at  some  optimum  pattern.  These  choices  are,  among  others: 

1.  20^  = 120®, 

2.  "triangular"  wave  envelope, 

3.  signal  length  = (T/2)  + €, 

4.  number  and  location  of  phase  changes,  and 

5.  position  of  signal  at  maximum  responses. 
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One  can  easily  see  the  many  combinations  that  must  be  studied  and  the 
effects  understood  if  a "good"  pulse  code  1b  to  be  found. 

Figure  5-10  displays  some  large  (-11  dB)  "grating"  lobes  at 
times  between  500  and  600  nsec.  An  example  of  the  pattern  tailoring 
that  can  be  accomplished  is  presented  in  Fig.  5-11.  In  that  pattern, 
the  late  "grating"  lobes  have  been  reduced  to  about  -19  dB  peak.  Hie 
coefficients  were  slightly  changed  in  two  respects:  (l)  the  time  of 

maximum  response  was  selected  to  be  earlier  than  in  Fig.  5-10,  and 
(2)  the  value  of  A^a)  was  changed  from  that  shown  in  Eq.  (5-47)  to 
a value  Ag(o()=sec  a.  It  can  be  seen  that  little  change  occurred  in  the 
major  lobe,  but  the  late  "grating"  lobes  were  reduced  8 dB. 
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VI.  IMPLEMENTATION- -SCANNED  AND  PREFORMED  BEAMS 


The  major  advantage  of  the  correlation  beamformer  is  its  ease  of 
implementation.  Delay  line  beamformers  are  limited  by  the  maximum 
length  of  the  delay  line  that  can  be  used.  Hence,  they  cannot  in 
general  take  advantage  of  the  diameter  of  the  cylindrical  array.  Phase 
shift  beamformers  have  not  taken  advantage  of  the  total  possible  aper- 
ture because  of  the  rise  time  problem.  Phase  shifters  also  require 
time  to  operate;  that  is,  several  cycles  of  the  signal  are  required  for 
an  analog  shifter  to  operate.  Compared  to  either  scanned  beam  operation 
or  preformed  beam  operation,  the  correlation  beamformer  can  scan  faster 
and  is  more  easily  implemented. 

The  following  sections  will  discuss  beam  scanning,  its  purpose, 
how  it  can  be  done,  and  at  what  speed.  Finally,  two  methods  of  imple- 
menting preformed  beams  will  be  presented. 

A.  Scanned  Beams 

The  ability  to  scan  a single  beam  at  a rate  which  has  no  theoretical 
upper  bound  is  perhaps  the  greatest  single  advantage  of  the  correlation 
beamformer.  The  scan  rate  is,  of  course,  limited  by  the  state  of  the 
art  in  electronics. 

Another  use  for  scanning  or  the  ability  to  rapidly  change  some 
parameter  of  the  formed  beam  is  in  time  varying  focus.  For  very  high 
resolution  arrays,  the  apertures  and  farfield  ranges  are  large.  If 
the  operating  range  is  substantially  less  than  the  farfield  range,  focusing 
at  the  expected  target  range  (a  function  of  time)  will  be  required.  The 
ability  to  rapidly  change  the  application  of  the  shading  coefficients  to 
element  output  voltages  is  a necessary  condition  to  time  varying  focus. 


93 


1. 


Method 


A repeat  of  Eq.  (3-1),  which  shows  the  essential  beamforming 
operation  with  the  output  carrier  intact,  indicates  the  factor  which 
leads  to  ease  of  implementation  and  to  the  unlimited  speed  of  scan. 

The  equation  is : 

N 

E*(cp,t)  = £ ARe  *(<p,t)  , (6-1) 

n=-N 


and  is  written  so  as  to  form  a beam  in  the  cpo=0  direction.  The  * 
indicates  a complex  wave,  or  in  real  life,  a real  wave  which  is  time 
varying.  Note  that  at  any  time  t'  infinitesimally  short  samples  of 
e^*  may  be  simultaneously  acquired.  These  samples  may  be  multiplied 
by  the  values  An  and  summed  to  derive  a sample  of  E*(q>, t').  With 
several  samples  of  E*,  the  beamformer  output  (including  carrier)  may 
be  reconstructed.  Note  that,  in  principle,  no  time  is  required  to 
sample  the  beam  output.  In  practice,  however,  the  sample  has  finite 
width,  and  the  multiply  and  sum  operations  require  time  as  determined 
by  the  bandwidth  of  the  circuits  performing  the  operation. 

Of  more  importance  in  terms  of  time  consuming  operations  is 
the  method  by  which  the  envelope  of  E*  is  determined.  Perhaps  the 
fastest  method  involves  quadrature  sampling  of  E*  and  retrieving 
the  envelope  via 

M = + (e2*)2  (6-2) 

where  the  samples  marked  1 and  2 are  taken  in  quadrature. 

Two  methods  are  available  for  determining  the  quadrature 
samples.  The  first  method  is  to  acquire  sample  E^*  at  any  time  t'  and 
wait  T/4  sec,  where  T is  the  carrier  period.  Then  E^*  is  acquired  at 
t'+T/4  in  the  same  way  the  first  sample  was  acquired. 
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by  \|c. 


The  second  method  is  to  use  two  sets  of  coefficients  differing 
Thus,  let 


2 

nA  = cos  a cos[2nR(l  - cos  a )] 

In  n n 

0A  = ccs  a cos[2jrR(l  - cos  a ) - \|c] 
2 n n n' 

and  then  perform  the  operations 


(6-3) 

(6-4) 


and 


Ej* 


N 

z 

n=-N 


,A  e * 
Inn 


i 


E2* 


N 

= z 

n- -N 


2Ane 


y 


where  e n*  is  only  one  sample  of  the  element  instantaneous  voltages 
taken  at  t',  but  ^ has  been  selected  to  shift  the  array  acoustic  center 
by  \/4.  In  this  way,  the  final  objective,  that  is,  a sample  of  |e*|, 
can  be  acquired  at  the  single  infinitesimal  time  t'. 


Now,  we  will  rewrite  Eq.  (6-l)  in  a form  more  descriptive 
of  the  scanning  operation.  Equation  (6-1)  is  written  so  as  to  form  a 
beam  in  the  cpo=0  direction.  The  coefficients  are 


In  general, 


p 

A (cp  =0)  = cos  (a  )cos[2itR(l  - 
n o ' n v 

to  form  a beam  in  the  cp  direction, 

o 

An(<P0)  = cos2(an-cpo)cos^2nR[l  - 


cos  a )] 
n 

we  have 
cos(an-cpo)]) 


(6-7) 


(6-8) 
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and  the  beamforming  equation  may  be  written  as 

N 

E*'<T-<P0,t)  = £ A (<PQ)  en*(<p,t)  . (6-9) 

n=-N 

Now,  in  order  to  scan  E*  as  a function  of  time,  we  simply  let 

cp0  = “s1  > (6-10) 

where  to  is  the  scan  rate  in  radians  per  second, 
s 

It  can  be  seen  in  Eq.  (6-9)  that  the  e*  samples  are  independent 
of  a-cp,  and  that  scanning  can  be  accomplished  by  simply  changing  the  An- 
The  fact  that  the  An  is  described  by  a continuous  function  of  a-<Po>  and 
that  the  function  is  even,  allows  the  continuous  function  A to  be  swept 
past  the  elements,  continuously  forming  a swept  beam.  Figure  6-1  shows 
a suggested  implementation  of  one-half  of  the  beamformer.  The  other  half 
contains  coefficients  of  the  ^A^  type  including  ♦,  so  that  by  combining 
their  outputs  using  Eq.  (6-2),  an  instantaneous  sweeping  estimate  of  the 
returned  energy  of  targets  can  be  obtained. 

2.  Scan  Speed  - Bandwidth  Requirements 


If  an  active  scanning  sonar  is  to  be  designed  with  BW  (degree) 
angular  beamwidth  and  range  resolution,  covering  the  azimuthal  plane 
out  to  a range  Pmax>  then  one  must  sample  56O/BW  beams  in  the  time 
required  for  the  transmitted  pulse  to  travel  Ap  yd.  Or,  since  range 
resolution  is  usually  given  in  t usee, 

t = fSe  , (6-11) 

ci 
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where  is  speed  of  sound  in  yd/ psec.  So,  one  must  scan  through  all 
360/BW  beams  in  t psec,  in  order  to  be  ready  to  repeat  the  operation 
during  the  next  range  increment  of  length  t.  Since  sonars  of  high 
resolution  (BW<0.5°)  may  contain  more  than  720  beams  and  have  a resolu- 
tion on  the  order  of  hundreds  of  microseconds,  it  is  most  important 
to  be  able  to  form  a beam  and  estimate  the  target  energy  in  a time  very 
short  compared  to  the  time  required  for  the  speed  of  sound  to  carry  one 
cycle  of  the  carrier  across  an  element.  To  our  knowledge,  the  correla- 
tion beamformer  is  the  only  process  capable  of  this  speed.  We  will  now 
develop  the  required  speed  or  bandwidth  for  the  multipliers  shown  by  X 
in  Fig.  6-1. 

A correlation  beamformer  with  total  sector  of  elements  of 
l80°  (2Q^=l80°)  can  form  a beam  with  width 

BW  = ^ , (6-12) 

where  K is  the  aperture  constant  in  degree-wavelengths . 

The  number  of  such  beams  needed  to  cover  the  azimuthal  plane  is 

= • (6-13) 

The  number  of  elements  required  to  cover  a cylinder  of  radius  R 
wavelengths  spaced  \/2  apart  is  ^nR,  2«R  of  which  are  used  in  any  one 
beam. 


Now,  if  720R/K  beams  are  to  be  scanned  in  t psec,  then  the 
coefficients  must  be  swept  past  ^«R  elements  in  r psec.  The  highest 
rate  of  change  of  coefficient  value  occurs  at  the  array  edge,  where 
the  rate  of  change  is  one  cycle  per  two  elements,  as  can  be  seen  from 
the  way  A^(a)  is  chosen.  Therefore,  the  highest  frequency  present  at 
the  coefficient  input  to  the  multipliers  is  found  to  be 
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FIGURE  6-1 

SCANNED  BEAM  MULTIPLIER  IMPLEMENTATION 
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Hz 


(6-14) 

Hz 

Of  course,  the  frequency  presented  to  the  element  input  of 
the  multiplier  is  just  f , the  transmitted  frequency,  a great  deal 
less  than  the  megacycle  rate  at  the  coefficient  input. 

Figure  6-2  shows  a plot  of  the  maximum  frequency  required 
for  the  multiplier  and  adder  in  a scanned  beam  sonar  of  beamwidth  EW 
degrees  and  range  resolution  equivalent  of  t psec. 


2nR 


rtK 

xBW 


In  order  to  test  the  simple  development  of  the  required 
frequency  response,  a fast  Fourier  transform  was  performed  on  a set 
of  shading  coefficients,  including  the  effects  of  finite  size  elements. 
The  results  are  plotted  in  Fig.  6-3,  and  show  the  time  function  at  the 
coefficient  input  to  the  multiplier  and  its  frequency  spectrum  through 
the  FFT.  The  time  scale  is  normalized  to  one  scan  per  second,  and  a 
small  radius  array  is  used  to  show  the  cyclic  nature  of  the  coeffi- 
cients. The  results  show  that  for  each  scan  per  second  at  R=12.5,  the 
maximum  frequency  at  the  multiplier  is  about  100  Hz.  For  one  scanner 
revolution  in  100  psec,  and  scaling  up  the  R from  12.5  to  300  (a  beam- 
width  of  0.1°)  the  maximum  frequency  is 


100 


300  1 sec 

12.5  ‘ 100  psec 


24.0  MHz 


(6-15) 


at  the  40  dB  down  point.  The  10  dB  down  point  in  the  spectrum  agrees 
very  closely  with  the  simple  development  given  earlier. 


Several  FFT's  were  run  for  various  R's  and  different  sectors 
of  active  elements,  with  and  without  finite  size  elements.  In  every 
case,  the  maximum  frequency  scales  directly  with  R.  Also,  the  case  of 


99 


BEAMWIDTH  - d.g 


) 


FIGURE  6-2 

RELATION  BETWEEN  SONAR  SYSTEM  RESOLUTION  PARAMETERS 
AND  THE  MAXIMUM  FREQUENCY  REQUIRED  FOR  SCANNER 
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finite  size  elements  (x=3-8j  see  Chapter  IV. A. 1.  and  Chapter  VIII. B.) 
exhibits  a flattened  spectrum  in-band  and  more  extended  tails,  with  the 
-20  dB  bandwidth  about  12^  wider  than  for  the  case  of  point  elements . 

A study  of  available  multipliers  indicates  that  frequency 
ranges  on  the  order  of  5 to  10  MHz  cam  easily  be  obtained,  and  20  MHz 
devices  are  available  at  reduced  but  adequate  accuracy.  This  would 
indicate  that  a scanning  sonar  can  be  constructed  to  operate  in  the 
0.3  to  0.1°  beamwidth  region  and  to  scan  the  entire  azimuthal  plane 
each  100  usee  (a  range  resolution  of  3 in.). 


B . Preformed  Beams 


If,  for  reasons  of  speed,  economics,  risk,  or  reduced  required 
sector  coverage  it  is  decided  to  commit  a separate  beamformer  to  each 
azimuthal  direction,  the  system  that  would  result  would  be  a preformed 
beam  system.  With  committed  beamformers,  peak  detection  and  post- 
detection filtering  can  be  accomplished  without  loss  of  time,  costing 
only  a time  delay.  This  section  will  outline  two  methods  of  implementing 
such  preformed  beams  excluding  the  envelope  detector. 

Again,  the  beamforming  operation  is  described  by  Eq.  (3-l)  or 
Eq.  (6-1) . The  coefficients  A^  are  scaling  factors  on  the  element 
voltages  e * which  are  to  be  sunned. 


1.  The  Resistor  - Operational  Amplifier  Scaling  Adder-Subtractor 

Figure  6-1*  shows  a circuit  of  the  resistor-operational  amplifier 
type.  When  is  adjusted  to  equalize  the  source  impedances  seen  by 
the  two  inputs  to  the  operational  amplifier,  the  dc  transfer  function 
for  the  circuit  is 


e 

o 


9 


(6-16) 
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where  Rp/R±i  is  the  magnitude  of  the  scaling  factor  for  e±i;  the  term 

is  added  for  e+^,  and  subtracted  for  e ^ . Equation  (6-l)  indicates 

that  terms  with  A X)  should  be  added  and  terms  with  A <0  should  be 
n n 

subtracted  from  the  total.  This  is  accomplished  if  the  terms  to  be 
added  correspond  to  e+1  and  the  terms  to  be  subtracted  correspond 
to  e . . 


Note  again  that  Eq.  (6-1)  forms  a beam  pointing  in  the  cpQ=0 

direction.  The  general  case  is  indicated  in  Eq.  (6-9),  which  can  be 

rewritten  in  a quantized  form,  indicating  the  output  of  the  jth  beam 

pointing  in  the  <p.  direction,  as 
J 

N+j 

E *(cp,t)  = £ A i(an-'Pi)  e *(<p,t)  , (6-17) 

J n= -N+j  n J n 

where  the  coefficients  are  computed  as  usual  using  the  argument  a -<p  , 

n 0 

and  an  points  to  the  nth  permanently  numbered  element  in  the  array. 

As  can  be  seen  by  the  argument  of  A,  if  the  element  spacing  and  beam 
spacing  are  equal,  then 


An,/an'V  = An+k,j+k^an+k"<PJ+k^ 


(6-18) 


for  any  k.  Also,  one  should  note  the  limits  on  the  sum  in  Eq.  (6-17). 
These  two  effects  result  in  the  symnetry  which  simplifies  the  implementa- 
tion of  many  preformed  beams.  That  is,  every  beamformer  is  identical 
to  every  other,  except  for  a unit  change  in  the  n for  the  contributing 
elements . 


If  the  element  spacing  is  smaller  than  the  beam  spacing,  some 
symnetry  still  exists.  For  example,  if  2Ap=3ACC,  then  two  types  of 
beams  can  be  formed,  and  alternately  placed  on  the  array. 
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A resistor-operational  amplifier  beamformer  containing 
97  resistors  has  been  constructed  and  used  in  the  TDA  tests  described 
in  the  next  chapter.  The  resistor-operational  adder  is  therefore  con- 
sidered low  risk.  The  system  does  require  the  use  of  preamplifiers  to 
drive  the  many  resistors  on  each  element. 

2.  The  Transformer  Adder 


Figure  6-5  shows  an  obvious  method  of  scaling  and  adding 
several  voltages.  If  all  transformer  primaries  have  equal  turns,  the 
scaling  can  be  accomplished  in  a quantized  manner  by  choosing  the  sense 
and  number  of  secondary  turns.  If  the  secondaries  are  in  series,  the 
voltages  add,  yielding  a quantized  version  of  Eq.  (6-9).  The  quantiza- 
tion occurs  on  the  A(n)  as  discussed  in  Chapter  IV. G.  and  proved  during 
the  experimental  validation,  Chapter  VIII. 

For  several  preformed  beams,  several  secondaries  can  be 
wound  on  the  same  transformer  core.  Some  experiments  were  conducted 
to  investigate  possible  problems  associated  with  the  transformer  adder. 
For  example,  it  was  shown  that  an  array  element  can  drive  a transformer 
primary  without  greatly  changing  the  element's  response.  Another  test 
involved  six  transformers,  each  with  20  secondaries.  Secondary  to 
secondary  feedover,  secondary  resonance,  and  null  voltages  were  measured. 
No  serious  defects  could  be  found  and  the  secondaries  combined  the 
transformer  outputs  just  as  expected.  Of  course,  the  more  transformers 
used  and  the  more  separate  secondaries  on  each,  the  more  feedover,  etc. 

It  was  not  possible  to  analyze  a large  matrix  of  transformers, 
except  in  the  simplest  way.  It  is  felt  that  a full  size  experimental 
beamformer  of  the  transformer  type  must  be  built  in  order  to  test  the 
feasibility. 
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FIGURE  6-  5 

TRANSFORMER  IMPLEMENTATION 
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The  advantage  of  the  transformer  beamformer  is  that  there 
are  no  active  devices  in  the  beamformer;  the  first  active  device  is 
the  beam  amplifier.  No  preamplifiers  mean  less  noise,  much  more  dynamic 
range,  and  higher  reliability. 


VII.  ANCILLARY  STUDIES 


During  the  course  of  the  study  of  the  correlation  beamformer, 
several  ancillary  studies  were  indicated.  This  chapter  describes  these 
studies  in  two  parts.  The  first  deals  with  the  practical  problems 
associated  with  array  construction.  The  second  part  will  describe  a 
short  but  useful  study  of  the  technique  of  linear  programming,  poten- 
tially useful  in  finding  codes  for  use  in  pulse  compression. 

A.  Array  Construction  Problems 

By  far  the  majority  of  work  reported  here  has  been  concerned  with 
beamforming  theory  and  computer  models.  In  every  case,  except  the  error 
analysis  (Chapter  IV.G.),  it  was  assumed  that  transducer  elements  generate 
a voltage  that  faithfully  reproduces  the  acoustic  pressure  on  the  element 
face.  The  elements  were  assumed  to  be  placed  in  a cylindrical  array 
with  perfect  accuracy. 

In  order  to  test  the  correlation  beamformer,  it  became  necessary  to 
construct  the  Two-Degree  Array  (TDA),  and,  to  test  the  pulse  compression 
codes  (Chapter  V),  it  will  be  necessary  to  construct  a much  larger  array. 
Several  practical  problems  had  to  be  addressed  in  order  to  build  high 
accuracy  arrays . 

1.  Element  Design 

Previously  constructed  arrays  at  ARL  have  consisted  of  staves 
of  several  elements  in  the  elevation  direction  (Z  direction  in  Fig.  4-2). 
The  purpose  was  twofold.  First,  when  the  dimension  of  any  side  of  a 
rectangular  parallelopiped  ceramic  element  becomes  much  larger  than  the 
resonance  dimension,  the  possibility  of  sympathetic  vibration  in  an 
unwanted  mode  exists . Such  unwanted  modes  can  adversely  affect 
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the  element  beam  patterns.  Second,  low  side  lobes  are  usually  required 
in  the  array  elevation  beam  pattern,  and  the  stave  elevation  pattern 
was  thought  to  be  the  main  contributor;  therefore  bizonal  shading  was 
applied  to  each  stave  of  several  elements  to  control  the  stave  elevation 
side  lobe  level. 

Our  discovery  that  a cylindrical  array  with  depth  in  the 
direction  of  signal  propagation  behaved  much  as  an  end-fire  array,  with 
sue  elevation  pattern  of  its  own,  led  to  the  attempt  to  simplify  the 
stave  design. 

If  a single  element  stave  (rectangular  parallelopiped)  i.° 
cut  to  be  resonant  in  the  thickness  mode  (the  intermediate  dimension), 
the  radiating  face  has  width  (small  dimension)  in  the  array  circum- 
ferential direction  and  length  ( long  dimension)  in  the  elevation 
direction.  The  width  is  controlled  by  the  required  element  spacing  on 
the  array  and  the  stave  length  is  an  open  parameter  which  can  be  used  to 
control  the  stave  elevation  response.  In  fact,  if  the  length  is  cut 
to  be  resonant  near  f /3,  the  stave  face  will  "flap",  superimposing  a 
sort  of  shading  on  the  face  velocity  function.  The  relative  amplitude 
of  the  long  mode  can  be  controlled  somewhat  by  mismatching  the  frequency 
relationship.  Figure  8-9  in  the  following  chapter  shows  an  element 
beam  pattern  from  the  TDA.  This  pattern  would  be  unacceptable  for  line 
arrays  because  of  the  large  second  lobes  at  -1^  dB.  For  the  cylindrical 
array,  however,  the  large  second  lobe  will  be  moderated  by  the  array 
pattern  (see  solid  curve  in  Fig.  3-7)-  The  first  lobe  will  be  nearly 
unaffected  by  the  array  pattern  because  of  its  nearness  to  the  major 
lobe.  The  composite  experimental  TDA  elevation  pattern  is  shown  in 
Fig.  8-10. 

By  taking  advantage  of  the  array  elevation  pattern,  it  was 
possible  to  construct  the  TDA  (next  chapter)  with  l8o  elements. 

Another  2°  beamwidth  array  constructed  earlier  at  this  laboratory 
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required  180  9-element  staves  ( 1620  elements  total),  but  it  was  no 
better  than  the  TDA  in  elevation  pattern  performance. 

2.  Element  Selection  Criteria 


During  any  array  construction  project,  elements  are  designed, 
the  design  is  tested  and  approved,  and  element  construction  begins.  In 
spite  of  all  precautions,  manufactured  elements  are  not  very  uniform 
in  performance.  Some  sort  of  element  selection  criteria  are  required 
to  choose  that  subset  of  elements  which  will  be  used  in  the  array.  The 
classic  criteria  have  been  a tight  grouping  of  frequency  at  which 
minimum  air  loaded  impedance  occurs . The  allowed  frequency  spread, 
and  in  fact,  the  very  choice  of  that  parameter,  appeared  to  be  more 
arbitrary  than  scientific. 

Actually,  one  would  hope  to  control  the  final  amplitude  and 
phase  response  of  each  element  after  mounting  in  the  array.  There  is 
some  question  that  "air  resonance  frequency"  is  a good  indicator  of  that 
amplitude  and  phase  control. 

During  our  construction  of  the  TDA  (an  experimental  array 
described  in  the  next  chapter),  elaborate  measurements  of  the  manu- 
factured elements  were  obtained.  Selection  of  the  elements  to  be  used 
in  the  array  was  accomplished  by  an  objective  comparison  of  air  impedance 
versus  frequency  plots  for  each  element.  After  construction  of  the  array, 
measurements  of  the  amplitude  and  phase  response  of  the  mounted  elements 
were  acquired.  A study  of  the  correlation  between  the  measurements 
made  in  air  before  mounting  and  measurements  made  in  the  array  after 
mounting  indicated  that  additional  in-air  measurements  (over  and  above 
"resonance"  frequency)  would  substantially  increase  the  uniformity  of 
elements  after  mounting. 

The  details  of  the  element  selection  criteria  and  the  method 
by  which  they  were  obtained  is  given  in  Ref.  9- 
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3-  Array  Construction  Technique 


An  additional  contributor  to  the  phase  scatter  of  mounted 
elements  in  an  array  is  the  inaccuracy  with  which  the  elements  are 
placed  on  the  cylinder.  The  normal  construction  technique  begins  with 
the  cylinder;  glue  and  backing  material  are  applied,  then  the  elements 
are  applied  to  the  backing.  Cylinder  inaccuracy,  glue  thickness  varia- 
tion, backing  material  thickness  variation,  and  element  thickness 
variation  all  combine  to  cause  the  element  faces  to  be  scattered  around 
the  proper  positions. 

A new  construction  technique  designed  to  eliminate  these 
construction  errors  has  been  developed  and  tested  on  a small  scale. 

The  results  indicate  that  the  method  can  increase  element  placement 
accuracy  by  a factor  of  ten  compared  to  the  old  methods . 

The  construction  method  involves  the  use  of  a mold,  which  is 
machined  to  hold  element  faces  in  the  proper  position.  Backing  material 
is  glued  to  the  element  backs  and,  finally,  the  cylindrical  baffle  and/or 
electronic  housing  is  inserted.  Hard  epoxy  is  used  to  connect  the 
element -backing  assembly  to  the  housing.  The  mold  is  then  removed, 
permitting  the  elements  to  be  wired.  Through  the  use  of  soft  epoxy 
and  a modification  to  the  mold,  the  face  material  is  poured. 

The  new  method  is  more  expensive  than  the  old  technique 
because  of  the  mold.  However,  the  extra  expense  can  be  justified  when 
large,  high  frequency  arrays  are  constructed  requiring  the  best  element 
placement  accuracy. 

B.  Linear  Programming  - An  Optimization  Technique 

Linear  programming  is  a mathematical  technique  by  which  an  objective 
function  can  be  optimized  while  constraining  other  functions  to  remain 


112 


within  bounds.  The  method  appears  to  be  useful  in  several  aspects  of 
the  correlation  beamformer.  For  example,  the  method  appears  to  be 
able  to  choose  or  solve  for  the  best  signal  code  for  use  in  the  pulse 
compression  problem.  Much  more  work  is  required  in  order  to  achieve 
that  capability,  however. 

Another  problem  or  question  involves  the  effectiveness  of  the 
heuristically  derived  shading  coefficient  selection  algorithms  described 
in  Chapter  III.  Can  other  coefficients  be  found  which  lead  to  better 
side  lobe  levels,  or  better  beamwidth?  How  nearly  optimum  are  the 
selection  algorithms?  In  order  to  answer  these  questions  and  learn 
more  about  the  technique  of  linear  programming,  it  was  decided  to  try 
to  linearize  the  beamformer  problem  and  exercise  the  supplied  software 
available  on  the  CDC  6600  computer  at  The  University  of  Texas  at  Austin 
campus . 

The  beamformer  problem  was  linearized  and  beam  pattern  runs  were 
compared  to  beam  patterns  from  the  algorithmic  approach.  The  results 
indicate  that,  for  R=12.3,  the  algorithm  approach  discussed  in  this 
report  is  nearly  optimum;  the  side  lobe  level  is  within  1 dB  of  the 
best  value  and  the  beamwidth  is  within  about  20$  of  the  best  value. 

The  optimum  coefficients  were  listed  and  plotted  but  no  sense  could  be 
made  of  the  computer  selection.  Certainly  no  scaling  rules  were  evident. 

It  appears  that  if  linear  programming  is  used  to  select  shading 
coefficients,  it  must  be  accomplished  for  the  specific  size  under  con- 
sideration. The  problem  cannot  be  run  for  very  large  arrays,  however, 
because  of  time  and  storage  requirements  of  the  software. 

Appendix  B gives  a brief  description  of  the  optimization  technique 
and  shows  the  beam  patterns  achieved  at  R=12.3.  It  is  felt  that  more 
work  on  this  subject  would  lead  to  a method  of  choosing  pulse  codes 
for  use  in  the  range  resolution  problem. 
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VIII . EXPERIMENTAL  VERIFICATION 


Two  separate  sets  of  experiments  have  been  conducted  during  the 
period  of  this  study  to  verify  the  predictive  capability  of  the  com- 
puter implemented  mathematical  model.  In  all  cases,  the  measurements 
confirmed  the  predictions,  not  only  in  basic  character,  but  also  in 
the  detail  structure  of  the  beam  pattern  side  lobes. 

The  first  measurements  were  conducted  on  an  existing  array  of 
very  small  radius  (iWO.  The  purpose  was  a "quick-look"  test  and  was 
performed  early  in  the  study.  The  second  tests  were  better  planned 
=md  utilized  an  array  especially  constructed  for  the  tests.  These 
experiments  will  be  described  separately. 

A.  AN/SQS-23  Array  Experiments 

In  early  April  197**  > the  STEP  barge  at  ARL's  Lake  Travis  Test 
Station  (LTTS)  became  available  for  use  on  this  contract.  An  AN/SQS-23 
sonar  array  had  been  suspended  under  the  barge  and  was  available  for 
testing.  The  array  is  physically  large  (8  ft  diam,  5 ft  height,  and 
5 ton  weight),  but  its  size  is  small  in  wavelengths  at  the  operating 
frequency  (R»M.  The  array  is  comprised  of  4^2  elements  (spaced  * 0.5  \), 
arranged  in  staves  of  nine  elements  each.  A contiguous  set  of  25  staves 
(the  front  half  of  the  array)  was  selected  and  routed  to  a test  beam- 
former  of  the  "variable  resistor -operational  amplifier  adder"  type. 

Because  of  the  large  physical  size,  the  experiment  was  rather 
limited,  and  because  of  the  small  normalized  size  (R*^),  the  beam 
patterns  do  not  look  "good".  However,  in  every  case  tested,  the  pre- 
dictions and  measurements  agreed. 
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The  test  included  a rough  error  sensitivity  measurement,  proved 
the  feasibility  of  steering  a beam  between  staves,  measured  both  sum 
and  difference  beam  patterns,  and  measured  sum  beam  patterns  as  a 
function  of  frequency.  In  every  case  tested,  the  comparison  between 
theory  and  experiment  was  very  good.  Figure  8-1  shows  an  example  plot 
of  one  of  the  measured  patterns . 

B.  The  Two-Degree  Array  (TDA)  Experiments 

During  October  1974,  tests  were  conducted  on  the  specially 
constructed  TDA  and  its  beamformer  test  bed.  The  system  had  been  under 
construction  since  early  in  the  contract,  and  its  purpose  was  to  make 
available  an  array  of  adequate  normalized  size  (R=l4.5)  for  detailed 
tests  of  the  method.  Of  course,  its  size  and  its  channel  bandwidth 
were  both  too  small  to  make  measurements  of  the  pulse  compression  ideas; 
however,  all  other  tests  were  to  be  performed  with  the  TDA.  Reference  10 
contains  a complete  description  of  the  system  and  the  tests  performed. 

In  this  section,  the  system  and  tests  will  be  sumnarized  in  order  to 
again  show  the  predictive  ability  of  the  computer  implemented  mathe- 
matical model. 

Figure  8-2  shows  the  completed  transducer  array.  It  contains 
180  hydrophones  spaced  on  2°  centers  on  the  cylinder  so  that  the  array 
outside  diameter  measures  17«53  in.  The  hydrophone  elements  are 
1.650  in.  by  0.624  in.  by  O.265  in.  blocks  of  CHANNELITE  5500® ceramic 
material,  resonant  near  100  kHz  in  the  0.624  dimension.  Single  element 
beam  patterns  were  measured  and  the  analytic  description  given  in 

Chapter  IV.A.l.,  Eq.  (4-5),  with  X=3.8  was  found  to  give  the  best  fit. 
Each  element  output  is  separately  conveyed  to  the  beamformer  test  bed, 
the  front  panel  of  which  is  shown  in  Fig.  8-3.  The  test  bed  panel  is 
composed  of  95*  10-turn  potentiometers,  whose  purpose  is  to  separately 
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control  the  shading  weight  of  the  array  elements  before  combining. 

In  this  way,  experiments  requiring  different  shading  coefficients  can 
be  easily  performed. 

The  experiments  were  conducted  in  the  ARL  sonar  model  tank,  which 

is  15  ft  wide,  60  ft  long,  and  12  ft  deep.  The  farfield  range  for  the 

2 

TDA,  estimated  from  D /\,  is  about  42  ft.  The  maximum  usable  test 

range,  however,  was  26  ft.  Thus,  the  majority  of  tests  were  conducted 

2 

at  0.62  D /\,  and  in  general  the  array  was  focused  for  that  range.  No 
true  farfield  beam  patterns  were  acquired. 

The  tests  performed  with  the  TDA  are  listed  and  the  results 
sumnarized  below. 

1.  Error  Sensitivity  Experiments.  Exact  shading  coefficients 
were  computed  and  linearly  quantized  into  several  levels,  ranging  from 
+1  to  -1  in  value.  Figure  8-4  shows  the  exact  coefficient  case,  and 
Fig.  8-5  shows  the  effect  of  quantizing  the  exact  coefficients  into 
five  levels  (+1,  +0.5*  0,  -0.5,  -l).  Both  cases  are  compared  to  the 
computer  model  for  that  case.  The  5-level  case  represents  the  most 
coarse  quantification  that  appears  useful. 

Other  interesting  patterns  were  obtained  in  this  category.  Examples 
include  a pattern  for  all  coefficients  uniform  (An=l)  and  the  case  of 
2-level  quantized  coefficients  An=+1  (see  Ref.  10)  . 

2.  Focusing  Experiments.  All  experiments  were  performed  at  some 

focused  range.  The  minimum  focused  range  used  was  12  ft.  This  range 

2 

represents  0.29  D /\  and  is  only  8.3  array  diameters.  Figure  8-6  is  a 
plot  of  the  pattern  for  a 12  ft  target  with  the  12  ft  focused  array. 

The  pattern  is  still  very  well  behaved.  The  pattern  resulting  from  a 
12  ft  target  when  the  array  is  not  focused  (or  farfield  focused)  is 
shown  in  Fig.  8-7.  As  expected,  without  focusing  the  sensitivity  is 
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down,  the  side  lobes  higher,  and  the  major  lobe  shape  is  bulbous. 

This  pattern  should  be  compared  to  Fig.  8-6. 

Many  other  patterns  were  acquired  to  determine  such  parameters  as 
depth-of-field. 

3-  Frequency  Response  Experiments.  The  array  was  matched  for  a 
particular  frequency,  then  several  patterns  were  acquired  at  different 
operating  frequencies.  The  purpose  was  to  determine  the  useful  bandwidth, 
and  if  the  computer  predictions  were  valid  for  this  case.  Figure  8-8 
is  presented  to  show  that  the  computer  prediction  is  good,  even  outside 
the  useful  bandwidth  of  the  correlating  processor. 

4.  Reduced  Sector  Experiments.  The  sector  of  active  elements 
was  reduced  to  8o°  and  then  to  40° . The  behavior  was  as  expected,  with 
computer  prediction  and  measurements  closely  correlated. 

5-  Split -Beam  Experiments.  Several  patterns  were  obtained  which 

represent  "half-beam"  and  "split-beam"  patterns.  The  patterns  presented 
herein  represent  the  sum  of  two  "half -beams".  The  "split-beam"  represents 
the  difference  of  two  "half-beams".  No  computer  predictions  were  made, 
but  the  behavior  was  as  one  would  intuitively  expect. 

6.  Elevation  Beam  Patterns.  Single  element  elevation  beam  patterns, 
of  which  Fig.  8-9  is  an  example,  and  the  array  elevation  pattern  through 
the  azimuthal  major  lobe  shown  in  Fig.  8-10,  were  obtained  and  compared 
with  the  theory  developed  in  Chapter  TV. A. 2.  As  can  be  seen,  the  rather 
unsightly  element  patterns  are  not  detrimental  to  the  array  elevation 
pattern;  the  latter  displays  side  lobes  more  than  25  dB  below  the  major  lobe. 

In  summary,  two  separate  experiments  were  performed  which  show 
that  the  method  is  a viable  beamforming  technique.  The  experiments 
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FIGURE  8-10 
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further  show  that  the  computer  implemented  mathematical  model  is  able 
to  predict  the  true  results,  not  only  in  basic  character,  but  also  in 
the  detail  structure  of  the  side  lobe  region. 
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IX.  CONCLUSIONS  AND  RECOMMENDATIONS 


It  has  been  shown  that  beam  formation  can  be  accomplished  with 
elements  or  staves  on  the  surface  of  a cylindrical  baffle  without  the 
use  of  delay  lines  or  pha§e  shift  networks.  The  shading  coefficients 
can  be  real  (vice  complex)  numbers.  The  ease  of  implementation  is  the 
greatest  advantage  accruing  from  the  use  of  the  method. 

This  report  has  summarized  the  results  of  a study  of  the  correlation 
beamformer.  Coefficient  selection  algorithms  have  been  developed  which, 
when  used,  result  in  a specific  useful  array  performance.  The  method 
may  involve  the  use  of  a coded  transmission  so  that  replica  correlation 
can  take  place  on  the  array  during  reception.  In  fact,  the  method  takes 
advantage  of  the  cylindrical  shape  of  the  array  in  all  respects  rather 
than,  as  was  done  before  this  study,  attempting  to  make  the  array  appear 
as  a plane. 

Another  advantage  of  this  approach  is  that  smaller  cylindrical 
arrays  can  be  made  to  deliver  equal  performance.  For  example,  the  TDA 
was  configured  to  deliver  equal  or  better  performance  in  a smaller  size 
than  the  cylindrical  transducer  (CT)  for  an  experimental  sonar.  The 
CT  was  2 6 in.  in  diameter  and  the  TDA  was  17  in.  in  diameter.  The  CT 
used  1620  elements  and  the  TDA  used  180.  Both  arrays  operated  at 
the  same  frequency  and  delivered  equal  beamwidths  with  the  same  side 
lobe  performance.  Appendix  C contains  other  examples  of  cost  and  size 
comparisons . 

The  method  is  relatively  easy  to  implement,  even  in  the  preformed 
beam  configuration.  Further,  Chapter  VI  describes  the  capability  of 
scanning  a single  beam  at  a rate  sufficient  to  cover  the  azimuthal  plane 
and  discern  targets  as  little  as  3 in.  apart.  Such  scanning  ability 
cannot  be  claimed  for  time  delay  or  phase  shift  beamformers.  Scanning 
a single  beam  further  reduces  the  hardware  required  to  construct  high 
resolution  sonar  soundheads. 
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Several  aspects  of  beamforraing  discussed  in  this  report  are  equally 
applicable  to  other  methods  of  beamforming.  For  example,  the  bandwidth 
and  the  elevation  beam  patterns  of  a phase  shift  beamforraer  should  behave 
with  array  depth  as  described  herein. 

It  is  recommended  that  the  following  steps  be  taken  to  further 
prove  the  method  and  to  show  conclusively  the  performance  and  the 
advantages  that  will  accrue  by  the  use  of  the  correlation  beamformer. 

1.  Only  one  code  class  for  pulse  compression  has  so  far  been  found. 
Although  it  does  accomplish  the  objective  of  pulse  compression,  it 
exhibits  several  undesirable  characteristics  which  should  be  eliminated. 
It  is  recommended  that  a test  of  the  presently  available  pulse  code  be 
accomplished  with  suitable  hardware.  The  test  should  include  angle  and 
time  response  measurements  of  the  array  and  tests  of  the  ability  of  the 
code  to  operate  in  bottom  reverberation. 

2.  It  is  recommended  that  an  additional  study  of  pulse  code6  be 
instituted  to  develop  a method  for  tailoring  the  array  time  response  to 
system  design  goals. 

5.  It  is  recommended  that  a full-scale  beam  scanner  be  designed  and 
constructed  to  verify  the  hardware  savings  that  can  be  achieved  with  the 
correlation  beamformer.  Such  a scanner  could  be  tested  with  the  array 
hardware  developed  in  Item  1. 

4.  It  is  recommended  that  suitable  hardware  be  developed  and  tested 
which  can  insonify  large  angular  sectors  with  the  pulse  code  transmission 
at  suitable  source  levels. 

5.  It  is  recommended  that  an  analytical  study  of  optimization 
techniques  (Chapter  VII  and  Appendix  B)  be  instituted.  The  technique  of 
linear  programming  has  already  demonstrated  its  capabilities  in  several 
beamformer  applications.  It  is  expected  to  be  especially  useful  in  pulse 
code  selection. 


152 


APPENDIX  A 

THE  ASSORTED  EFFECTS  OF  REDUCING  THE 
SECTOR  OF  ACTIVE  ELEMENTS 


135 


APPENDIX  A 


THE  ASSORTED  EFFECTS  OF  REDUCING  THE 
SECTOR  OF  ACTIVE  ELEMENTS 

The  reference  array  considered  in  Chapter  III  of  this  report  contains 
active  elements  on  one-half  of  the  cylinder  (20^180°).  The  array  depth, 
in  the  direction  of  signal  propagation,  is  therefore  R (the  array  radius 
in  wavelengths).  This  array  depth  has  been  shown  to  influence  or  control 
the  elevation  beam  patterns,  the  frequency  response  of  the  beamformer, 
and  the  array  rise  time. 

For  whatever  reason,  if  the  sector  of  active  elements  is  reduced, 
there  is  an  attendant  reduction  in  the  array  depth,  the  array  width 
(the  aperture),  the  number  of  contributing  elements,  and  the  maximum 
dimension  of  the  array.  If  some  parameter,  such  as  beamwidth,  is  to 
remain  unaffected  as  the  sector  reduction  takes  place,  the  problem 
becomes  more  complicated. 

This  appendix  will  present  design  equations  and  plots  of  their 
functional  dependence  which  will  be  useful  in  studying  or  designing  a 
reduced  sector  array. 

A single  beam  requires  an  aperture  W on  a cylinder  of  radius  R, 
which  results  in  an  array  depth  L.  The  aperture  Is  composed  of  a sector 
of  active  elements  covering  an  angle  on  the  cylinder  of  20^.  The  design 
may  include  several  preformed  beams  covering  a region  in  space  measured 
by  the  angle  S.  These  statements  are  illustrated  in  Fig.  A-l.  As  can  be 
seen  from  Fig.  A-l, 

W = 2R  sin  (A-l) 

and 

L = R(1  - cos  ctjj)  . (A-2) 
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If  one  desires  a beamwidth  of  BW  (degrees),  then  an  aperture  width 
of  W=k/BW  wavelengths  is  required,  and  R must  grow  to  keep  pace  with 
any  decrease  in  sin  O'  . That  is, 


R = 


R F, 
o 1 


(A-3) 


where 


and 


R 


o 


W K 

2 ' 2 BW(deg) 


(A-4) 


F,  = 


1 sin 


“n 


(A-5) 


a function  which  is  shown  graphically  in  Fig.  A-2. 

If  one  wishes  an  array  depth  of  L wavelengths,  and  the  array  radius 
is  R,  then  either 


L = R(l  - cos  ctjj) 


(A-6) 


or 


n L 

(l  - COS  Ojj) 

(A-7) 

Both 

of  these  functions  are  plotted  in  Fig.  A-3. 

One  can  compute  the  area  of  the  resolution  cell  of  an 
reduced  sector  by  noting  that 

array  with 

£y  = pBW(rad)  = ^ ^So 

(A-8) 

and 

AX  = ~ ~ (1  - cos  dtjj)  , 

(A-9) 

137 


Function 


FUNCTIONS  - F<  AND  F 


or 


R \ 

AX  = - — 7 — — (l  - cos  a ) 
2 sin  a.  ' N 

N 


(A- 10) 


if  the  beamwidth  is  held  constant. 


By  normalizing  to  the  resolution  (AXq)  which  is  obtainable  when 


2a  =l80° , we  have 
N 


R X. 

- Hr  ’ 


(A-ll) 


so  that 


1 - cos 


AX  = AX 


V 


3 in 


= HF2 


(A -12) 


This  function  also  describes  the  relative  reduction  in  resolution  area, 
1 • e • , 


1 - cos 


V 


A = Ao\  ■ sin  aN  7 ■ AoF2 


(A- 13) 


since  Zfy  is  unaffected.  The  function 


1 - co.  “h 

F2  " .1.0,  * tan  T 


is  plotted  in  Fig.  A- 4. 

Another  consequence  of  reducing  the  sector  of  elements  to  form  a 
single  beam  is  the  possible  increase  in  the  total  required  hardware 
(element  count)  for  an  array  of  several  preformed  beams. 

If  the  angular  sector  width  of  active  elements  for  one  beam  is  2a^ 
radians,  and  if  the  sector  width  covered  by  additional  preformed  beams 
is  S radians,  then  the  angle  subtended  by  elements  on  the  cylinder  is 
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20jj+S  radians.  The  required  length  of  arc  on  the  cylinder  is  R(2ctN+S); 
there  are  two  elements  per  wavelength,  so  the  number  of  elements, 
preamplifiers,  and  filters  each  is 

M = 2R(2aN+S)  , (A- 15) 


or  in  general 


2R 

o 

sin  Ojj 


(2aN+s) 


(A- 16) 


If  we  divide  the  general  part  count  by  the  standard  part  count,  we  have 

a number  F^  that  relates  the  fractional  increase  in  hardware  caused  by 

the  reduction  of  the  active  sector  2a„:  thus, 

N ’ 


M(aN,S)  (2aN+S) 

= Mo(i’s)  = u+s)  sin  “N 


(A-17) 


This  function  is  plotted  in  Fig.  A-5-  One  can  see  that  for  only  one 
preformed  beam  (S=0),  a hardware  savings  results  as  2a^  is  decreased. 
However,  as  the  sector  covered  by  preformed  beams  increases,  there  exists 
an  cijj  below  which  one  cannot  work  without  significantly  increasing  the 
required  hardware  - 


An  additional  function  F^  relates  the  largest  dimension  LD  of  the 
array  to  changes  in  2a^.  The  sector  of  active  elements  is  2otjj+S.  The 
chord  of  this  arc  is 


LD  - 2R  sl"(“B  * |)  for  aB  * | « | 
LD  = 2R  otherwise. 


(A-18) 


Again,  to  preserve  BW,  replace  R by  RQ/sin  a^,  so  that 


2R 

o 

sin  Qtjy 


lU2 


(A-19) 
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By  normalizing  to  the  2Q^=l80°  case,  we  have 


LD 

2R 


sin 


K + D 


sin  a, 


N 


for  + | 5 | , (A-20) 


or 


LD  _ 1 

2R  ~ sin 


aN 


= F, 


for 


“»  + f >f 


(A-21) 


A plot  of  F^  is  presented  in  Fig.  A-6. 


In  order  to  summarize  and  outline  the  use  of  these  functions, 
consider  the  following  design  example: 

BW  = 0.3°  , 

f = 300  kHz  , 

S = 120°  = y rad  (400  beams),  and 
t = pulse  length  = 100  psec, 


resulting  in  a theoretical  range  resolution  of 

AX  = ~ = P000  x 1(^°  x 10.  . ■ = o.25  ft 


First,  the  array  is  designed  with  2aN=l80°,  and  the  array  parameters 
are  evaluated  as  follows: 


(1) 

\-t0~  3^7500  " 0’01667  n 

(2) 

Bo  - W - ■ 98-53  • 

R X = 98.33  X 0.01667  = 1.64  ft 

* 
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(3)  Mo  = 2R0(n+S)  = 2 x 98-33  x (it  + y)  = 1030  , 

(4)  AX„  . ^ g O.O1667  , 0.81958  ft  , 

(5)  LDo  = 2Rq  = 2 X 98.33  = 196.67  , 

or 

ldq\  = 196.67  x 0.01667  = 3.278  ft 

From  line  (4)  above  we  have  the  range  resolution  of  O.819  ft,  and  we 
wish  to  achieve  a resolution  of  0.25  ft.  Since 

AX  = AXqF2 

and 

F2  = Ol9  = °'5°5 

it  can  be  computed  from  Fig.  A-k  that  2aN  must  be  68°  to  achieve  the 
desired  resolution. 

We  are  now  prepared  to  modify  each  of  the  parameters  found  in  lines 
(2)  through  (5)  through  the  use  of  20^=68° : 

(6)  R = R F. 

o l 

With  reference  to  Fig.  A-2,  compute 

f1(68°)  = 1.788  , 

so  that 

R = 98.33  x 1.788  = 175.8  , 

or 


R\  = 175.8  x 0.01667  - 2.93  ft 


I 


M = M F, 
o 3 

With  reference  to  Fig.  A-5>  compute 

F^So^S)  = F5(68°,  120°)  = 1.12 

so  that 

M = 1030  X 1.12  = 1155 

(8)  LD  = 2R  Fj, 

o 4 

With  reference  to  Fig.  A-6,  compute 

F4(68°,  120°)  =1.788  , 

so  that 

LD  = 2 X 98.33  X 1.788  = 351.68  , 

or 

LD\  = 351.68  X O.OI667  = 5.86  ft 

Thus,  the  adding-only  array  for  these  parameters  will  have  the  following 
design  features: 

maximum  dimension  = 5.86  ft  , 

number  staves  = 1155  , 

(preamplifiers,  etc.) 

range  resolution  = 0.25  ft  , and 

(equivalent  to  t=100  pisec) 

radius  of  curvature  of  array  = 2.93  ft 

It  should  be  noted  that,  even  though  reducing  the  sector  20^  will 
solve  some  range  resolution  problems,  the  best  solution  is  to  find  a 
suitable  pulse  shape  that  will  give  the  range  resolution  by  pulse  com- 
pression. When  compression  is  accomplished,  this  array  can  be  built  in 
a 3.3  ft  package,  compared  to  the  5.9  ft  indicated  above. 
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BEAMPORMER  OPTIMIZATION  TECHNIQUES 


A.  Introduction 


At  the  present  time  there  exists  no  direct  analytic  techniques  for 
the  selection  of  shading  coefficients  to  form  beams  with  elements  placed 
in  a nonplanar  arrangement.  Design  techniques  in  use  today  are  either 
approximations  to  the  planar  case,  some  algorithmic  method  as  described 
in  this  report,  or  else  some  iterative  technique.  The  technique  of 
"phasing-to-a-plane"  is  an  example  of  approximating  a plane  array. 

Trial  and  error  is  an  example  of  an  iterative  technique.  A somewhat 
more  sophisticated  iterative  technique  is  linear  programming,  which  is 
the  principal  subject  of  this  appendix. 

Linear  programming  has  been  applied  to  the  problem  of  determining 
the  real  coefficients  for  a circular  array  of  given  element  number  and 
spacing.  The  coefficients  were  determined  by  this  method  to  obtain  as 
closely  as  possible  a desired  array  beam  pattern. 

The  results  of  this  linear  programming  investigation  have  been 
encouraging  from  several  aspects.  First  of  all,  the  validity  of  the 
design  technique  discussed  in  other  sections  of  this  report  has  been 
confirmed.  For  example,  those  design  procedures  are  seen  to  produce 
beam  patterns  within  1 dB  in  peak  side  lobe  response  and  within  about 
20^  of  the  minimum  beamwidth  attainable  with  the  linear  programming 
optimum  design.  Second,  however,  this  investigation  has  shown  that  there 
are  some  worthwhile  gains  to  be  made  by  the  application  of  optimization 
techniques  to  beamformer  design.  Besides  the  narrower  beamwidth  and 
improved  side  lobe  control,  there  are  indications  that  the  response  of 
the  array  to  a pulse  signal  can  be  tailored,  within  limits,  to  meet  other 
system  design  goals. 
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The  study  reported  here  was  brief  and  centered  on  only  one 
optimization  technique,  linear  programming.  There  are  other  optimization 
methods  which  are  perhaps  more  suited  to  the  problem  at  hand.  Linear 
programming  is  not,  strictly  speaking,  applicable  to  the  beamformer 
optimization  problem.  Therefore  we  have  approximated  the  true  problem 
by  another  problem  which  is  directly  solved  with  the  linear  program. 

Some  other  practical  optimization  technique  may  be  more  directly  appli- 
cable to  the  beamformer  optimization  problem.  This  should  be  the 
subject  of  a further  investigation. 

The  following  sections  detail  the  linear  programming  approach  to 
beamformer  optimization. 

B.  Linear  Programming 

The  linear  programming  problem  that  we  will  solve  can  be  stated  in 
the  following  form:  find  the  {a^},  i=0,l,2, . . . ,N,  subject  to  the 

constraints 

N 

S A c..  S b , j = 1,  2,  ...,  M (B-l) 

i_0  1 J1  J 

M £ N + 1 

such  that 

N 

£ A.  d (B-2) 

i=0 

is  maximized.  We  do  not  require  the  A^O  as  usually  stated  in  the  more 
general  linear  programming  problem. 

One  characteristic  of  linear  programs  is  that,  if  a solution  exists, 
it  can  be  found  in  a finite  number  of  iterations.  There  are  also  pro- 
cedures for  determining  if  a solution  exists,  or  if  the  solution  is 
unbounded. 
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The  primary  reason  for  using  linear  programming  to  optimize  the 
fceamformer  is  that  there  are  many  computer  programs  available  to  solve 
the  problem.  This  limits  the  beamformer  optimization  problem,  so  far  as 
we  are  concerned,  to  setting  up  the  problem  in  the  form  previously  stated 
and  then  using  the  appropriate  routine  to  find  the  array  coefficients. 

In  the  next  section  we  discuss  the  approximation  of  the  true  beamformer 
problem  by  one  that  fits  the  desired  form. 


C.  The  Optimization  Problem 


The  output  of  the  circular  array  to  a signal  incident  from  direction 
cp  can  be  expressed  as  the  convolution  of  the  array  impulse  response  and 
the  input  signal: 


h(t,cp) 


(B-3) 


where 

A = coefficients  of  the  2N+1  elements, 

I(an-cp)  = individual  element  beam  pattern, 
s(t)  = signal, 

c = sound  speed,  and 

a^,  cp,  and  r are  as  defined  by  Fig.  3-1*  The  time  reference  for  both  s(t) 
and  h(t,cp)  are  taken  to  be  the  center  of  the  array.  Here  h(t,cp)  is  the 
response  of  the  array  to  the  signal  s(t)  incident  on  the  array  in  the  cp 
direction,  and  so  h(t,cp)  actually  represents  a directivity  function  or 
beam  pattern  for  the  array. 


If  the  A are  constrained  to  be  symmetric  so  that  A =A  , n=l, 
n n -n 

2,...,N,  and  the  individual  element  beam  patterns  are  also  symmetric 
about  their  facing  directions,  then  the  resulting  beam  patterns  will  be 
symmetric  about  cp=0  and  we  have 


h(t,cp)  = AQl(cp)s(t  + ^ cos  cp) 


N ( 

2 I(V 

n=l  { 


cp)s[t  + £ cos(an-cp)] 


+ l(an+cp)s[t  + £ cos  (a 


n-«p)]|  • 


(B-4) 


Since  the  time  of  arrival  of  the  signal  is  ordinarily  unknown,  a 
more  useful  characterization  of  the  array  beam  pattern  is  the  envelope 
of  h(t,cp),  which  we  designate  by  H(t,cp). 


H(t,cp) 


-E 


h2(t,cp) 


(B-5) 


where  the  positive  square  root  is  taken  so  H(t,cp)  is  nonnegative  and 
h ' (t,cp)  is  the  Hilbert  transform  of  h(t,cp)  with  respect  to  the  time 
variable  t.  It  is  easy  to  show  that 


h’(t,cp)  = AQI(cp) s ' (t  + ^ cos  cp) 


+ 


I(an-cp)s'[t  + ^ cos(an-cp)] 


(B— 6) 


+ I(a  +cp)s'[t  + - cos(a  -Kp) ] 
n c n 

where  s'(t)  represents  the  Hilbert  transform  of  s(t).  We  note  that  both 
h(t,cp)  and  h'(t,cp)  are  linear  functions  of  the  |Anf,  but  H(t,<p)  is  not. 

Now,  in  order  to  form  a beam  pattern  with  peak  side  lobe  level  SLL 
at  time  t=tQ,  we  require 


H(tQ,Cp)  =1,  <p  = 0 , 

(B-7) 

H(to,cp)  S SLL,  ®b  S cp  < 100°  , 

(B-8) 

w 
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where  cp^  is  the  angle  at  which  the  main  lobe  response  first  reaches  the 
peak  side  lobe  level.  Thus  2cp^  is  the  total  width  of  the  main  lobe  at 
the  side  lobe  level.  The  response  at  <p=0 , the  direction  of  the  main 
lobe,  is  required  to  be  unity. 


Another  approach  to  the  requirement  stated  by  Eq.  (B-8)  is  to 
quantize  the  angle  cp  in  the  range  (cp^»l80°).  If  the  quantization  interval 
on  cp  is  small  enough,  then  Eq.  (B-8)  will  effectively  be  met.  Then  the 
side  lobe  constraints  become 


H(tQ,<p  ) 5 SLL  J = 1,  2,  ...,  M 

and  the  main  lobe  response  can  be  assured  by  maximizing 

H(tQ,q>),  9-0 


(B-9) 


(B-10) 


Equations  (B-9)  and  (B-10)  then  are  similar  to  Eqs.  (B-l)  and  (B-2) 
except  that  Eqs.  (B-9)  and  (B-10)  are  not  linear  functions  of  the  {Anf. 
However,  the  constraints  expressed  by  Eq.  (B-9)  will  be  met  if  all  the 
following  linear  constraints  are  met  simultaneously: 


h(t  ,<p  ) 5 0.707  SLL 

® J 

-h(to,q>  ) § 0.707  SLL 
h'(t  ,<p  ) 5 0.707  SLL 

J 

-h'(to,«Pj)  S 0.707  SLL 


) J = 1,  2, 


J 


• • • } 


M 


(B-11) 


This  set  of  Eq.  (B-ll)  overly  constrains  the  beam  pattern,  since  Eq.  (B-ll) 
is  a sufficient  but  not  necessary  condition  for  Eq.  (B-9).  By  a similar 
line  of  reasoning,  Eq.  (B-10)  can  be  replaced  by  maximizing 

h(t  ,cp)  9 = 0 (B-12) 


We  have  thus  linearized  the  problem  by  substituting  the  set  of  4m 
constraints  Eq.  (B-ll)  for  the  M constraints  Eq.  (B-9),  and  by  changing 


the  objective  function  Eq.  (B-10)  to  Eq.  (B-12).  Equations  (B-ll)  and 
(B-12)  are  in  the  form  specified  for  Eqs.  (B-l)  and  (B-2);  therefore 
the  linear  programming  routine  may  be  applied  to  the  solution  of 
Eqs.  (B-ll)  and  (B-12). 

In  the  next  section  some  of  the  results  of  the  linear  programming 
approach  to  the  problem  using  Eqs.  (B-ll)  and  (B-12)  are  presented. 

D.  Linear  Programming  Results  and  Conclusions 

The  linear  programming  technique  discussed  in  the  previous  section 
was  applied  to  the  design  of  a 77-element  array  of  elements  spaced  on 
half-wavelength  centers  around  a semicircle  of  radius  12.5  wavelengths. 
The  objective  of  the  design  was  to  obtain  a main  lobe  beamwidth  as  narrow 
as  possible  while  maintaining  uniformly  low  side  lobes. 

Several  slightly  different  parameter  values  were  chosen  for  test. 

The  parameter  values  used  and  the  most  important  design  values  resulting 
therefrom  are  given  in  Table  B-l.  The  parameter  cp^  controls  the  main 
lobe  width  since  the  first  side  lobe  constraint  is  placed  at  that  angle. 
The  remaining  side  lobe  constraints  were  placed  at  1°  intervals  from  5° 
to  90° • The  parameter  tQ  is  used  to  control  the  phasing  of  the  signal 
to  which  the  array  is  "matched"  by  the  optimization  process.  In  all 
these  designs  the  signal  was  a sine  wave  with  constant  amplitude.  The 
value  to=-0.05  corresponds  to  the  situation  where  a peak  of  the  sine  wave 
signal  is  coincident  with  the  front  of  the  array.  For  tQ=0.0,  the  signal 
peak  was  l8°  behind  (toward  the  array  center)  the  front  of  the  array. 

The  latter  phasing  gives  slightly  better  results,  but  no  effort  was  made 
to  determine  the  optimum  phasing. 

For  reference,  the  beam  pa -tern  obtained  by  choosing  the  element 
gains  according  to  the  rules  outlined  in  this  report  is  given  in  Fig.  B-l, 
and  the  design  values  are  also  given  in  Table  B-l.  The  coefficients 
chosen  were  for  the  reference  array. 


156 


TABLE  B-l 


COMPARISON  OF  BEAM  PATTERN  PARAMETERS  BETWEEN  ALGORITHMIC 
APPROACH  AND  LINEAR  PROGRAMMING  TRIALS 


TEST 

DESIGN 

PARAMETERS 

PEAK 

SIDE  LOBE 
(dB) 

BEAMWIDTH  (deg) 

^b 

t 

o 

-3  dB 

At  Peak. 
Side  Lobe 
Level 

0* 

— 

— 

-16.2 

2.5 

5.3 

1 

2.7 

-0.05 

-17.0 

2.2 

4.5 

2 

2.0 

-0.05 

-15.1 

2.0 

3.9 

3 

2.7 

0 

-16.7 

2.2 

4.5 

k 

2.2 

0 

-16  A 

2.1 

4.3 

* 

Test  0 is  the  reference  array  case . 

(Test  0 and  1 are  shown  in  Figures  B-l  and  B-2,  respectively.) 
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FIGURE  B - 1 

AZIMUTHAL  BEAM  PATTERN 
USING  REFERENCE  ARRAY  COEFFICIENTS 


ARL  - UT 
AS-  76-  56 
DAS  - DR 
1-27-76 


In  all  of  the  designs  shown  in  Table  B-l,  the  individual  element 
patterns  were  assumed  to  be  [(l+cos  g)/2r  , where  the  angle  p for  each 

element  is  measured  with  respect  to  the  direction  that  element  is  facing. 
Better  beam  patterns  not  included  in  this  report  were  obtained  assuming 
unity  gain,  semicircular,  individual  element  patterns. 

In  the  design  values  measured  and  listed  in  Table  B-l,  it  is 
apparent  that  the  linear  program  designs  are  somewhat  better  than  the 
design  of  Fig.  B-l.  Main  lobe  widths  approximately  14$  narrower  at  both 
the  3 dB  and  the  peak  side  lobe  level  points  have  been  obtained  without 
giving  up  peak  side  lobe  level,  thereby  improving  the  directivity  index 
by  almost  1 dB.  Perhaps  more  important,  the  troublesome  first  side  lobe 
in  Fig.  B-l  has  been  effectively  suppressed. 

It  should  be  pointed  out  that  the  beam  patterns  achieved  in  the  four 
examples  in  the  table  are  almost  certainly  not  optimal  because  the  side 
lobes  are  of  varying  height.  Figure  B-2  shows  an  example  pattern  from  the 
linear  program.  A better  design  would  reduce  the  largest  side  lobes  at 
the  expense  of  raising  the  smaller  side  lobes  to  yield  an  optimal  Chebyshev 
pattern  in  which  all  the  side  lobes  are  the  same  height.  Such  a design 
can  be  approached  with  the  repeated  application  of  the  techniques  dis- 
cussed here  but  no  successful  optimum  designs  were  accomplished  during 
this  short  study.  4 

There  are  other  important  beamformer  parameters,  bandwidth  and 
sensitivity,  for  instance,  which  have  not  been  directly  investigated  in 
this  brief  study.  It  does  not  appear  that  linear  program  designs  will 
prove  to  be  difficult  because  of  these  considerations.  In  fact,  the 
linear  programming  approach  provides  a means  of  directly  addressing  many 
such  problems. 

There  remains  one  area  of  difficulty  in  the  direct  application  of 
linear  programming  to  beamforraing  as  discussed  here.  That  is  in  dealing 
with  large  arrays  of  many  elements.  While,  in  principle,  the  technique 
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is  independent  of  size,  practical  considerations  limit  the  size  of  the 
arrays  that  may  be  handled  to  about  40  wavelengths  radius  and  250 
elements.  Perhaps  larger  arrays  could  be  dealt  with  on  the  very  large 
computers  available  today  but  the  cost  would  be  high.  The  patterns 
reported  here  for  a 12.5  wavelength,  77-element  array  were  generated  in 
about  4 min  each  on  a CDC  6600.  The  computer  size  and/or  time  require- 
ments grow  exponentially  with  problem  size. 

The  solution  to  this  difficulty  lies  not  in  larger  and  faster 
computers,  but  in  a better  algorithm  to  perform  the  optimizations.  The 
linear  program  algorithm  is  probably  not  a very  good  technique  for 
solving  the  beamformer  problem;  however  it  is  readily  available  and 
therefore  quite  useful  in  an  exploratory  program.  A really  serious 
attempt  at  the  optimization  of  large  arrays  of  the  type  considered  in 
this  study  will  require  a more  fundamental  look  at  optimization  methods. 
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APPENDIX  C 

A SIZE  AND  COST  COMPARISON 

During  the  past  several  years  ARL  has  been  involved  with  either 
the  specification,  design,  or  construction  of  at  least  three  small, 
high  resolution  sonars,  whose  application  and  requirements  represent 
the  solution  to  real,  practical  problems.  Because  of  the  applicability 
of  these  systems,  they  were  selected  for  use  in  a realistic  comparison 
with  the  new  method. 

The  correlation  method  derives  its  advantage  from  the  single 
fact  that  no  time  delay  (delay  lines)  or  phase  shift  networks  are 
required.  The  beamformer  cost  is  therefore  reduced.  In  addition, 
because  no  practical  delay  line  constraints  are  imposed,  the  entire 
diameter  of  a cylindrical  array  may  be  used  as  the  aperture,  resulting 
in  either  smaller  arrays  for  the  same  beamwidth  or  higher  resolution 
in  the  same  sized  arrays. 

One  should  keep  in  mind  that  the  three  systems  used  for  comparison 
are  already  compromises  that  Eire  favorable  to  the  delay  line  technique 
and  size  constraint.  It  is  therefore  fair  to  ask  such  questions  as 
how  the  methods  would  compare  if  the  angular  resolution  was  doubled  or 
the  sector  covered  by  preformed  beams  increased. 
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The  first  system  used  for  comparison  is  called  the  Cylindrical 
Transducer  (CT).  It  represents  the  experimental  application  of  a 
cylindrical  array  to  an  existing  sonar  system.  The  array  and 
beamformer  simultaneously  form  60  2°  beams  to  cover  120°  in  azimuth, 
and  are  electrically  switchable  in  120°  increments  for  a full  360° 
overall  coverage.  The  total  array  diameter  is  26  in.  The  beamformers 
use  60  delay  lines  costing  about  $70  each,  for  a total  of  about  $4.2K. 
Each  azimuthal  beam  is  formed  using  50  contiguous  staves,  for  a total 
of  60  beams;  this  requires  5000  shading  resistors  costing  about  $300. 

The  application  of  the  correlation  beamformer  to  this  system 
would  use  all  the  staves  on  the  front  half  of  the  array  for  the 
aperture,  thereby  reducing  the  required  array  diameter  to  17  in. 

However,  9°  of  the  total  180  staves  would  contribute  to  each  beam, 
increasing  the  number  of  shading  resistors  to  90x60=5400,  costing  about 
$540.  No  delay  lines  would  be  required,  nor  would  additional  staves 
or  preamplifiers/filters . Therefore,  the  cost  savings  would  be 
$4200- $240=$ 3960*  The  array  size  would  be  reduced  to  about  65$  of  the 
comparable  CT  array. 

The  second  system  used  for  comparison  is  called  the  VSS  HAWAII. 

The  array  which  is  cylindrical  and  19*5  in.  in  diameter,  forms 
forty- five  2°  beams  covering  a 90°  sector.  There  is  no  requirement 
for  additional  sector  coverage,  so  that  only  86  staves/preamplifiers/ 
filters  are  required.  The  VSS  was  implemented  with  45  delay  lines  at 
about  $50  each  (about  $2.25K).  Since  42  staves  contribute  to  each 
of  45  beams,  the  system  uses  I890  shading  resistors  (about  $189). 

If  a correlation  beamformer  uses  a l80°  sector  of  active  elements 
that  are  contributing  to  each  beam  (an  array  size  of  11.5  in.),  and 
45  beams  are  required,  the  total  number  of  staves/preamplifiers/filters 
would  have  to  be  increased  to  126.  The  VSS  uses  86.  The  increased  cost 
of  preamplifiers/filters  in  the  correlation  implementation  at  about 
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|15  each  would  be  $600.  Also,  with  84  staves  contributing  to  each  of 
45  beams,  3780  resistors  will  be  required,  which  cost  more  than  the  VSS 
implementation  by  about  $190.  However,  the  savings  by  not  using  delay 
lines  is  $2.25K.  Therefore,  the  correlation  beamformer  could  perhaps 
be  constructed  for  about  $1.46k  less  than  the  delay  line  implementation. 
Note  also  that  the  size  is  reduced  59$. 

This  system  represents  a difficult  comparison  because  of  the 
restricted  requirements  of  sector  coverage.  If  the  requirement 
was  to  cover  360°  with  2°  beams,  the  same  number  of  staves/preamplifiers/ 
filters  would  exist  in  both  implementations.  The  differences  would  then 
be  a saving  of  I9K  in  delay  lines  and  a factor  of  two  increase  in  the 
number  of  resistors  costing  about  $756,  for  a net  saving  of  $8.244K. 

These  cost  savings  depend  almost  linearly  upon  the  required  azimuthal 
coverage,  with  the  savings  being  greatest  at  full  coverage.  On  the 
other  hand,  the  delay  line  method  is  cheapest  if  only  one  beam  is 
formed . 

The  third  comparison  concerns  the  AN/VQS-1.  This  sonar  array  is 
cylindrical  with  a diameter  of  14  in.  and  contains  54  staves  covering 
approximately  170°  on  the  cylinder.  The  AN/^QS-1  forms  40  beams  each 
3°  wide  covering  a 120°  azimuthal  sector.  The  system  contains  54 
preamplfiers  costing  about  $13  each  ($702).  Only  15  staves  contribute 
to  each  of  40  beams,  requiring  600  resistors  ($60)  and  40  delay  lines 
at  $50  each  ($2K).  The  total  cost  attributable  to  the  beamforming 
method  is  $2.J6k. 

The  correlation  method  applied  to  this  beam  requirement  can  be 
accomplished  by  placing  100  staves  over  300°  of  a cylinder  whose 
diameter  can  be  5.8  in.  compared  to  14  in.  for  the  AN/^QS-1.  More 
preamplifiers  and  resistors,  which  will  be  required,  will  cost  $1.54K. 
This  represents  a size  savings  of  41$  and  a cost  savings  of  44$  ($1.22K). 
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Another  way  to  compare  the  methods  is  to  hold  the  size  (l4  in.), 
operating  frequency,  and  cost  ($2.76k)  constant,  and  determine  what 
resolution  can  be  achieved  over  which  sector  of  beam  coverage  by 
using  the  correlation  method.  The  results  of  such  computations  show 
that  bearawidths  of  1.24°  can  be  achieved  on  a 14  in.  cylinder.  Beams 
may  be  added  to  cover  a 40°  sector  (32  beams)  before  the  cost  exceeds 
the  |2.76k  figure.  If  120°  coverage  is  required,  then  the  cost  would 
be  f4.>48K;  for  full  360°  coverage  using  1.24°  beams,  the  cost  would 
be  I7-975K. 

In  summary,  the  correlation  method  is  generally  cheaper  to  apply 
and  results  in  a smaller  cylindrical  array  when  compared  to  other 
methods.  Depending  primarily  upon  the  sector  coverage,  the  cost 
savings  may  be  significant  or  insignificant;  however,  in  all  cases 
tested,  the  size  reduction  was  nearly  5 0 
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APPENDIX  D 

SIGNAL- TO- NOISE  RATIO  ESTIMATES 

During  a progress  meeting  in  April  1975  concerning  the  comparative 
advantages  of  the  correlation  beamforraer  versus  a phase  shift  beamformer 
(delay  line  implemented),  the  subject  of  processing  gain  or  signal  to 
thermal  noise  ratio  arose.  A few  quick  calculations  indicated  that  the 
S/N  of  the  correlation  beamformer  would  be  6 to  12  dB  less  than  the 
phase  shift  implementation. 

Upon  returning  from  the  meeting,  a relatively  quick  analysis  was 
conducted  which  was  released  as  a technical  letter.  Both  the  cover 
letter  (because  it  contains  the  summary)  and  the  technical  letter  are 
reproduced  here  since  they  represent  the  total  effort  on  S/N  analysis 
conducted  during  this  study. 

Both  the  correlation  method  (called  Dave's  method)  and  the  phase 
shift  beamformer  (called  Bobby's  method)  were  being  compared  for 
application  to  a specific  requirement.  The  broadside  line  array, 
inserted  as  a standard,  did  not  itself  meet  the  system  requirements, 
but  served  as  a convenient  reference. 
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Enel:  (l)  "Signal  to  Noise  Ratio  or  Processing  Gain  Against 

Thermal  Noise  for  Different  Beamforming  Techniques," 
TL-EA-75-2  of  29  April  1975 


Our  conversation  last  Thursday  concerning  the  relative 
performance  of  Bobby's  phase  shift  beamformer  and  my  replica- 
correlation  beamformer  has  preyed  on  me  since  that  time.  As  I 
recall,  your  opinion  was  that  my  method  would  result  in  Signal 
to  Thermal  Noise  Performance  between  6 and  12  dB  poorer  than 
Bobby's  phase  shift  method. 

Your  argument,  simple  and  elegant  as  it  was,  convinced  me 
also,  and  I found  myself  agreeing  on  Friday,  "certainly  not  12  dB, 
but  perhaps  6 to  ^ dB  loss!" 

I have  computed  both  signal  to  thermal  noise  ratios  since  my 
return  home  and  find  that,  in  fact,  the  difference  between  the  two 
methods  is  0_.1*5  dB  to  my  detriment.  Enclosure  (l)  should  enable 
you  to  come  to  the  same  conclusion. 

I admit  to  some  haste  in  this  draft;  however,  several  of  my 
crew  have  checked  these  figures,  and  they  appear  to  be  in  order. 

I'm  sure  you  can  understand  my  interest  in  correcting,  as  soon  as 
possible,  the  assertion  at  the  Friday  meeting  that  my  method  was 
defficient  in  signal  to  noise  ratio. 
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By  the  way,  when  compared  to  the  broadside  line  array  with 
equal  aperture  and  equal  element  weights,  it  is  true  that  my  method 
(and  Bobby's)  compare  unfavorably;  however,  the  broadside  line  does 
not  meet  the  specification.  An  equivalent  line  would  contain  1*00 
elements,  and  as  shown  in  Section  I of  enclosure  (l),  the  reference 
numbers  are : 


Signal  Gain  = 1*00 
Noise  Gain  = 20 

Signal/Noise  Processing  Gain  = 20  +26  dB 

Section  I of  enclosure  (l)  outlines  the  ground  rules  and 
assumptions  made  in  this  analysis.  It  also  solves  the  broadside 
line  reference  case  for  comparison  purposes.  Section  II  details 
the  analysis  for  Bobby's  method,  and  Section  III  computes  the  same 
number  for  my  method.  The  results  are  given  in  the  following  table. 
In  each  case  shown,  the  apertures  are  equivalent. 


NUMBER  OF  ELEMENTS  PER  BEAM 

ARRAY  SIGNAL 
Gain 

Gain  (dB) 

ARRAY  N0IS§ 

Gain 

Gain  (dB) 

S/N  RATIO 
Number 
dB 


REFERENCE 


BOBBY 

DAVE 

BROADSIDE 

256 

60o 

1*00 

190.92 

+U5. 62 

153.89 

+U3.71* 

Uoo.oo 

+52.014 

12.76 

+22.12 

10.83 

+20.69 

20.00 

+26.02 

lU.96 

+23.50 

lU.20 

+23.05 

20.00 

+26.02 

As  you  can  see,  my  results  indicate  that  altho’igh  my  method 
is  less  sensitive  and  exhibits  less  S/N  than  the  reference  case  of 
the  broadside  line  (which  does  not  meet  specs),  both  Bobby  and  I 
have  the  same  problem.  It  seems  to  be  the  price  paid  with  either 
method  for  meeting  the  required  specifications.  Please  note  the 
insignificant  difference  between  the  S/N  ratios. 

I appreciate  your  interest  and  criticism. 


David  A.  Smith 
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SIGNAL  TO  NOISE  RATIO  OR  PROCESSING  GAIN  AGAINST 
THERMAL  NOISE  FOR  DIFFERENT  BEAMFORMING  TECHNIQUES 


David  A.  Smith 


I.  ARRAY  GAIN 

SIGNAL  TO  THERMAL  NOISE  RATIO 

The  purpose  of  this  section  is  to  identify  and  describe  the  ground  rules 
for  computation  and  comparison  of  the  Array  Signal  Gain,  the  Array  Thermal 
Noise  Gain,  and  the  Processor  Signal  to  Thermal  Noise  Ratio  for  the  two 
different  processors  under  consideration  and  a standard  broadside  line  array. 

In  one  case  the  processor  is  a phase  shift  beamformer,  which  I will 
call  Bobby's  method.  In  the  other  case,  the  replica-correlation  beamformer, 
which  I will  call  Dave's  method. 

I will  assume  that  each  element  in  the  arra>  has  identical  receiving 
sensitivity  and  an  identical  thermal  noise  source.  Further,  for  the  present 
case  I assume  omnidii  ectional  element  patterns  so  that  each  ..ement  generates 
the  same  peak  signal  output  voltage  e when  insonified  with  a standard  sound 

O 

field.  The  thermal  noise  voltage  generated  by  each  element  is  similarly 
e^  (an  RMS  noise). 

I also  assume  that  no  additional  noise  is  introduced  by  the  processor 
itself,  and  that  any  channel  bandwidth  factor  lowering  the  noise  is  assumed 
to  be  incorporated  in  e . 

Now,  the  Signal  Gain  is  a number  which  relates  the  processor  output 
signal  (peak)  S to  the  peak  eleir  nt  voltage  eg.  It  may  be  stated  as  either 
a numerical  factor  or  in  dB  (20  log). 
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The  Thermal  Noise  Gain  is  similarly  the  ratio  of  beamformer  output  noise 

N to  the  element  output  noise  e„ . Both  numbers  are  in  RMS  volts  and  are  due 
o N 

only  to  thermal  sources  in  each  element. 

The  Signal  to  Noise  Ratio  is  the  ratio  of  peak  signal  to  RMS  thermal 
noise  at  the  output  of  the  processor  to  that  ratio  produced  at  the  output 
of  each  element.  That  is,  it  is  the  processor's  signal  to  noise  ratio  gain 
or  Processing  Gain  against  thermal  noise. 

In  both  processors,  there  exist  a set  of  shading  coefficients  A^  which 
are  used  to  scale  the  element  output  before  summing  or  before  phase  shifting 
and  summing. 

Now,  in  general,  the  SIGNALS  combine  coherently;  that  is,  we  know  the 

"bh 

phase  of  each  element's  signal.  Let  the  i — element  deliver  e^(t)  volts 

where  e.(t)  is  the  instantaneous  output  voltage  when  the  peak  is  e„.  Long 
1 b 

duration,  co-sinusoidal  signals  are  assumed  so  that 

e.  (t)  = e cos[cp(t,i)]  (l-l) 

1 b 

If  each  element  voltage  e^(t)  is  scaled  by  A^  (the  shading  coefficient) 
and  then  summed  over  the  M contributing  elements,  we  have  the  processor  output 
signal  S . That  is , 

M 

S = A^e^(t).  (l-2) 

i=l 

In  Bobby's  case,  the  A^  contain  phase  shifts  which  correct  for  the  elements 
position  in  the  array.  In  Dave's  case,  the  A^  are  real  numbers  that  represent 
the  signal  replica.  These  separate  analyses  appear  in  Sections  II  and  III. 

The  Thermal  Noise  sources,  each  generating  e^  (RMS  volts),  are  combined 
in  a Root  Mean  Square  sense  after  scaling  by  A^ , because  the  individual  phases 
are  unknown  and  random.  Therefore,  the  processor  output  noise  N^  is 
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N = 

o 


M 


E (v»)! 


i-1 


1/2 


(1-3) 


Again,  separate  analysis  of  the  two  methods  appear  in  Sections  II  and  III. 


The  element  output  signal  to  noise  ratio  or  processor  input  signal  to 

noise  ratio  is  — . The  processor  output  signal  to  noise  is  rr-  . The 
*"N  N 

Processing  Gain  is  therefore  S_  fN  ° 

b O 

An  Example 

The  standard  reference  case  is  the  broadside  line  of  M equal  spaced 

elements  of  equal  strength.  Such  a line  would  exhibit  a far field  beam 

pattern  like  x . 

x 

Fran  Eq.  1-2,  the  signal  is 


S = 


M 

22  Ai6i(t) 


(1-U) 


For  the  broadside  line  with  target  in  the  farfield,  all  elements  receive 
the  same  signal,  and  if  t is  selected  so  that  the  peak  signal  is  measured, 
and  equal  strength  shading  means  Ai  = 1 for  all  i,  we  have 


(1-5) 


(1-6) 


From  Eq^  1-3,  the  noise  is 
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N = 
o 


M 


£ (W 


i=l 


1/2 


(1-7) 


Again  A = 1,  so  that 


M 


i=l 


2 2 
leN  = eNM 


(1-8) 


or  N , 

-s-  = yr 


(1-9) 


The  Processing  Gain  against  thermal  noise,  for  this  reference  case, 
is  therefore 


(1-10) 
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II.  BOBBY'S  METHOD 


The  Phase-Shift  Processor 


From  Section  I,  Eq.  1-2,  we  have  that  the  processor  output  signal  S 
is 

M 

S = Ai  6i(t)  (2_1) 

i=l 

where  the  A?  contains  the  phase  shift  information.  Since  from  Eq.  1-1, 
ei (t ) = eg  cos  [cp(t,i)] , 

M 

S = } * A*  eg  cos  [ 9(  t , i ) ] (2-2) 

i=l 

But  the  complex  part  of  A?  is  designed  to  he  just  enough  to  bring  each 
element  output  voltage  in  phase,  so  that 


M 

S = A.  eg  cos  (cp(t)) 

i=l 


(2-3) 


(where  the  A^  are  now  real  numbers),  and  if  t is  chosen  to  be  the  time  of 
cos  (<p(t))  = 1,  then  we  have 


( 2-U ) 


(2-5) 
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Data  Analysis 


The  values  of  Bobby's  shading  coefficients  are  functionally  very 
complicated,  and  have  been  computed  and  the  output  listed  on  the  CDC  3200 
computer.  The  data  extracted  axe  noted  "best"  as  of  2 6 April  1975. 

Table  2-1  indicates  the  heading  information,  and  Table  2-2  lists  the  real 
part  of  the  shading  coefficients.  One-half  are  listed,  since  the  array 
is  even  symmetric. 


TABLE  2-1 

HEADING  INFORMATION  FOR  "BEST"  COEFFICIENTS 


Number  of  Elements 
Frequency 
Array  Radius 
Stave  Height 
Stave  Width 
Stave  Spacing 
Total  Angular  Width 
Velocity  of  Sound 


256 

300  kHz 
36  inches 
1 in. 

0.150  in. 

0.010  in. 

6U.9^  deg. 
58,267.6  in/sec 
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TABLE  2-2 


BEST  ALS  COEFFICIENT  AS 


1 

2 

3 

.99998 

.99998 

•99957 

.99793 

•99711 

.99630 

.99262 

.99100 

.98937 

.98369 

.98166 

.97925 

.97161 

.96880 

.96600 

.956U6 

.95291 

. 91*936 

.93800 

.931*09 

.92981+ 

.91675 

.91217 

.90759 

. 89287 

.88766 

.8821+5 

.86625 

.86075 

.851*97 

.83739 

.83136 

.82532 

.80650 

.80000 

•79351 

.77360 

.76690 

.76001 

.73917 

.73211 

.72505 

.7031*3 

.69609 

.68875 

.66651 

.65906 

.65152 

.62882 

.62121 

.61360 

. 59060 

.58289 

•57519 

.55206 

.5M<3l* 

.53663 

.51351 

.50582 

.1+9812 

.1*7516 

.1*6756 

.1*5996 

.1*3731 

. 1+2979 

.1*2235 

.1*0011* 

.39280 

.3851+7 

.36381 

.35673 

.31*965 

. 32868 

.32171* 

. 311*96 

1.01625 

1.73111 

2.1*1*598 

OF  26  APRIL  75 


1*  5 


•99916 

.99875 

5 

•99507 

.99385 

10 

.98771* 

.98572 

15 

• 97683 

.971*1*1 

20 

.96282 

.95961* 

25 

.91*581 

.9l<190 

30 

.92559 

.92131* 

35 

.90268 

.89778 

1*0 

.8772U 

.87171* 

1*5 

.81*920 

.81+31*2 

50 

.81905 

.81278 

55 

.78701. 

.78031 

60 

.75312 

.71*623 

65 

.71781* 

.71063 

70 

.6811+1 

.67396 

75 

.61*398 

.6361*1* 

80 

.60593 

.59826 

85 

.5671*9 

•55978 

90 

.52892 

.52120 

95 

.1*901+7 

.1*8281 

100 

.1*5236 

.1*1*1+81* 

105 

.1*11*91 

.1*071*7 

110 

.37825 

.37103 

115 

.31*256 

.33562 

120 

.30817 

.30139 

125 

128 

J.62 


A program  was  written  for  the  HP-65,  the  purpose  of  which  was  to  calculate 

g 

the  values  represented  in  Eq.  2-6  for  — , and  Eq.  2-9  squared  for  the  value 

es 


GT 


The  results  of  the  computation  with  the  special  program  was: 


— = 95.H60U2 

eS 


(2 


— = 81. U 0017 


(2 


M = 128 


(2 


Because  256  elements  will  be  used  and  the  array  is  even  symmetric, 
the  values  were  multiplied  by  2 and  the  square  root  of/N  was  obtained 


— = 190.92  (2 

eS 

N 

— = 12.76  (2 

eN 

M = 256  (2 

In  order  to  check  the  computations,  a standard  program  for  the  HP-65 
was  used  (Mean,  Standard  Deviation,  Standard  Error  - STD-02A).  Equation  2-6 
can  be  rewritten  in  terms  of  the  average  A.  = A as 

f—  - MA  . (2. 

es 

Further,  Eq.  2-9  can  be  recast  in  terms  of  the  average  A^  = A and 
the  standard  deviation  of  A^  = *K  as  follows: 
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jL 


-10) 

-11) 

-12) 

-13) 

■11*) 

•15) 

■16) 


1 


N / 2 _2 

— =\l  *A  (M-l ) + M A 
eN 


S No 

Thus  by  computing  A and  , one  can  evaluate  — and  — . 

eS  eN 

The  values  resulting  from  this  check  computation  are  as  follows: 
A = O.T>+578 
= 0.28350 
M = 256 

Therefore,  from  Eq.  2-l6,  we  have 

— = MA  = 190.92  , 

es 

and  by  Eq.  2-17,  we  have 
N 

— = 12.76  . 
eN 

These  values  will  be  seen  to  agree  with  the  previous  results. 

In  summary,  I present 


Array  Signal  Gain  (Bobby's  Method) 
— = 190.92  4 +1*5.62  dB 


Array  Noise  Gain  (Bobby's  Method) 
N 

— = 12.76  4 +22.12  dB 

eN 

Processing  Gain  (Bobby's  Method) 

— =£  = 1U.96  4 +23.50  dB 
6S  ”o 
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(2-17) 


(2-18) 

(2-19) 

(2-20) 

(2-21) 

(2-22) 


(2-23) 

( 2-2U ) 
(2-25) 


III.  DAVE'S  METHOD 


The  Replica-Correlator  Processor 

From  Section  I,  Eq.  1-2,  we  have  that  the  processor  output  signal  S is 
M 

S=X  Ai6i(t) 

i=l 

where  the  A.  are  real  numbers  and  no  phase  shift  is  implied.  Since  from 
Eq.  1-1, 

ei(t)  = es  cos  [cp(t,  i ) ], 

the  value  of  S becomes 
M 

S =2  Vs  COS  1 ) 3- 

i=l 

Now,  in  Dave's  Method,  the  A^  can  be  written  as 

A = cos  [2jtR(1-cos  a.)]  cos2  a 
1 i i 

where tJ  is  the  array  radius  in  wave  lengths  and  a.  is  the  pointing  angle  to 
the  i element.  If  the  elements  are  spaced  Q wavelengths  apart  on  the  cyl- 
inder. 


So 


Ai  = 


(3-1) 


(3-2) 


(3-3) 


(3-4) 


cos  [2*R(l-cos(i§))]cos2(!§) 
R R 


(3-5) 


Now  the  signal  eg  cos  (cp(i,t))  is  of  the  same  form  as  part  of  the  A.; 
that  is,  if  t is  selected  so  that 

es  cos  [cp(i=0,t=to)l  =eg  (3-6) 

then  the  signal  at  the  i—  element  is 

6i  = 6S  C0S  C2irR(1_cos^))]  (3-7) 

If  we  substitute  Eqs.  3-5  and  3-7  into  Eq.  3-3,  we  have  that 
M 

eS  cos2  C2irR  ( l-c°s  (■^^) ) ]c°s2  (-^^)  • (3-8) 

i=l 


Let 

Al.  = cos  [2itR  (1-cos  (I2))! 

R 

and 


A 2. 

1 


,iQx 
cos(— ) 


> 


and  rewrite  Eq.  3-8  as 
M 


i T,  (“iV 


i=l 


(3-9) 


(3-10) 


(3-11) 


Now  for  the  noise  factor,  from  Eq. 


1-3  we  have 


(3-12) 


or 


(3-13) 
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Data  Analysis 


Since  simple  functional  forms  exist  for  A1  and  A£L,  a FOCAL  program 
was  written  for  the  FDP-8  to  evaluate  Eq.  3-11  and  the  square  root  of  Eq.  3-14. 
The  program  listing  is  presented  in  Table  3-1.  -• 

For  the  case  of  601  elements  spaced  0.5236  \ apart  on  l80°of  a 100  \ 
radius  cylinderical  array,  the  results  are  as  follows: 

Array  Signal  Gain  (Dave’s  Method) 

- = 153.89  =»  +43-74  dB 

es 


Array  Noise  Gain  (Dave's  Method) 
N 

— = 10.83  =»  +20.69  dB 

SN 


Processing  Gain  (Dave's  Method) 

S e 

— r— - = 14.21  =>  +23.05  dB 

eq 

O O 
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TABIE  3-1 


C-FOCAL,  1971 

01.05  C ADDISON  SIG.  4 NOISE  FACTORS 

01.10  S P=3.l4l59 

01.40  A R,N,S, !!! 

01.70  S Y=2*P*R 

02.10  S SG=0;S  N0=0 

02.30  F I=-N,N;S  X=I*S/R;D  3;D  4 

02.40  S NS=FSQT(NO) 

02.50  S SN=SG/NS 

02 . 60  T 2*N+1, SG, NS, SN, ! ! ! ;Q 

03.10  S A1=FC0S ( Y* ( 1-FCOS ( X ) ) ) 

03.20  S A2=FC0S(X) 

04.10  S SG=SG+A1*A1*A2*A2 

04.20  S N0=N0+(Al*A2*A2)t2 
* 


INPUTS  WERE  R=100  N=300  (M  IN  TEXT  = 2N+1)  ELEMENT 
SPACING  S(Q  IN  TEXT)  = 0 .5236. 


references 


1.  Walter  Welkcwitz,  "Directional  Circular  Arrays  of  Point  Sources," 

J.  Acoust.  Soc.  Am.  28,  362-366  (1956). 

2.  T.  T.  Taylor,  "Design  of  Line-Source  Antennas  for  Narrow  Beam- 
width  and  Low  Side  Lobes,"  IRE  Trans . on  Antennas  and  Propagation, 
pp.  16-28,  January  1955* 

3.  C.  L.  Dolph,  "A  Current  Distribution  for  Broadside  Arrays  Which 
Optimizes  the  Relationship  Between  Beam  Width  and  Side-Lobe  Level," 
Proc.  of  the  IRE  on  Waves  and  Electrons,  pp.  335-3 48,  June  1946. 

4.  R.  C.  Spencer,  "Fourier  Integral  Methods  of  Pattern  Analysis," 

Report  762-1,  Radiation  Laboratory,  Massachusetts  Institute  of 
Technology,  Cambridge,  Massachusetts  (January  1946). 

5.  Henry  J.  Bickel,  "Real-Time  Spectrum  Analysis,"  J.  Sound  and  Vib., 
14-20  (March  1971)- 

6.  C.  L.  Temes,  "Sidelobe  Suppression  in  a Range-Channel  Pulse- 
Compression  Radar,"  IRE  Trans,  on  Military  Electronics,  pp.  162-169, 
April  1962. 

7.  D.  G.  Tucker  and  B.  K.  Gazey,  Applied  Underwater  Acoustics  (Pergamon 
Press,  New  York,  19 66),  Chapter  6. 

8.  Bernard  D.  Steinberg,  "Comparison  Between  the  Peak  Sidelobe  of 
the  Random  Array  and  Algorithmically  Designed  Aperiodic  Arrays," 
teee  Trans,  on  Antennas  and  Propagation,  pp.  366-370,  May  1973* 

9.  Michael  Kump,  David  A.  Smith,  and  Robert  M.  Adams,  "Multiple  Factor 
Selection  Technique  for  Ceramic  Transducer  Elements,"  Applied 
Research  Laboratories  Technical  Report  No.  75-68  (ARL-TR-75-68), 
Applied  Research  Laboratories.  The  University  of  Texas  at  Austin, 
Austin,  Texas  (6  January  1976). 

10.  Bruce  0.  Moses  and  David  A.  Smith,  "Acoustic  Receive  Beamforming 
Experiments  Using  a Cylindrical  Array  with  Real  Shading  Coefficients," 
Applied  Research  Laboratories  Technical  Report  No.  74-56  (AHL-TR- 
74-56),  Applied  Research  Laboratories,  The  University  of  Texas  at 
Austin,  Austin,  Texas  (19  December  1974). 


191 


October  1976 


DISTRIBUTION  LIST  FOR 
ARL-TR-76-2  (Revised) 

FINAL  REPORT  UNDER  CONTRACT  N61339-74-C-0031,  Task  0003 
29  November  1973  - 31  December  1975 
UNCLASSIFIED 


Copy  No. 

Commanding  Officer 
Naval  Coastal  Systems  Laboratory 
Department  of  the  Navy 
Panama  City,  PL  32401 


1 - 

3 

Attn:  Code  112.2 

4 - 

13 

Code  112.6 

14 

A.  L.  Rolle 

15 

J.  L.  Bealor 

16 

Don  Folds 

Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 


17 

Attn: 

C.  E.  Fox 

18 

R.  Staith,  Code  663CB 

19 

C.  D.  Sknlth,  Code  06HI 

20 

W.  L.  Thompson 

21 

A.  Breece,  Code  661C-12 

22 

R.  Snuggs,  Code  660F-5 

23 

M.  T.  Giver,  Code  660E-23 

24 

V.  Graves,  Code  06H4-3 

25 

W.  Brandon,  Code  06H4-42 

26 

LCDR  D.  F.  Bolka,  Code  06H2 

27 

J.  E.  Neely,  Code  06H1-1 

28 

M.  C.  Hansen,  Code  06HB 

29 

E.  J.  McKinney,  Code  0342 

30 

F.  J.  Romano,  Code  03423 

31 

B.  Welch,  Code  03424 

32 

Roy  Heaney,  Code  66lT 

33 

D.  H.  Bennet,  Code  FMS-4o6B 

Commander 

Naval 

Undersea  Center 

Department  of  the  Navy 

San  Diego,  CA  92132 

34 

Attn: 

J.  R.  Reardon,  Code  NUC- 35601 

35 

H.  W.  Volberg 

36 

R.  A.  McLennan,  Code  NUC  256 

37 

D.  L.  Carson,  NUC- 302 

38 

Mike  Kolar,  NUC  6512 

193 


V 

Distribution 

List  for  ARL-TR-76-2  (Cont'd) 

Copy  No. 

39 

40 

41 

Officer- in-Charge 

New  London  Laboratory 

Naval  Underwater  Systems  Center 

Department  of  the  Navy 

New  London,  CT  06320 

Attn:  H.  Newman 

W.  R.  Schumacher 
W.  Konrad 

42 

Commanding  Officer 

Newport  Laboratory 

Naval  Underwater  Systems  Center 

Department  of  the  Navy 

Newport,  RI  02840 

Attn:  T.  M.  Fitzgerald 

43 

44 

Director 

Naval  Research  Laboratory 
Department  of  the  Navy 
Washington,  DC  20375 
Attn:  John  Munson 

L.  Van  Buren 

45 

46 

Commanding  Officer 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Deptvrtment  of  the  Navy 
Silver  Spring,  MD  20910 
Attn:  Bill  Rust 

Joe  Mobley 

47 

Commanding  Officer 

Naval  Facilities  Engineering  Command 

Department  of  the  Navy 

200  Stovall  Street 

Alexandria,  VA  22332 

Attn:  P.  Cave,  Code  NAVFAC-032C 

48 

Commanding  Officer 

Naval  Air  Systems  Command 

Department  of  the  Navy 

Washington,  DC  20361 

Attn:  L.  E.  Wheat,  NAVAIR-5305A 

I 


* 


Distribution  List  for  ARL-TR-76-2  (Revised)  (Cont'd) 
Copy  No. 

49  Office  of  the  Chief  of  Naval  Material 

Department  of  the  Navy 
Washington,  DC  20360 

Attn:  CDR  R.  K.  Hastie,  NAVMAT  0314 

50  Superintendent 

Naval  Postgraduate  School 

Monterey,  CA  93940 

Attn:  Librarian  (Code  2124) 

51  Commanding  Officer 

Naval  Electronic  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 
Attn:  W.  Kamminga,  NAVELEX-3203 

52  Commander 

Naval  Ship  Engineering  Center 
Department  of  the  Navy 
Center  Building,  Prince  George's  Center 
Ryattsville,  MD  20782 
Attn:  F.  Hill,  NAVSEC  6126 

Chief  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA  22217 

53  Attn:  R.  Obrochta,  Code  ONR-222 

54  M.  Blizard,  ONR-212 

55  Commanding  Officer 

Naval  Air  Development  Center 
Department  of  the  Navy 
Warminster,  PA  18974 
Attn : Bill  Parlglan 

56  Naval  Research  Laboratory 
Underwater  Sound  Reference  Division 
P.0.  Box  8337 

Orlando,  FL  32806 
Attn:  Joe  Blue 


195 


■ ■■ 


Distribution  List  for  ARL-TR-76-2  (Rev.)  (Cont'd) 


Copy  No. 

57 


58 


59  - 70 


71 


72 

75 

74 

75 

76 

77 

78 

79 

80 

81 

82 


Catholic  University 
Electrical  Engineering  Department 
620  Michigan  Avenue,  NE 
Washington,  DC  20017 
Attn:  Dr.  H.  F.  Harmuth 

Defense  Advanced  Research  Projects  Agency 
Tactical  Technology  Office 
lkOO  Wilson  Blvd. 

Arlington,  VA  22209 
Attn:  CAPT  H.  Cox 

Commanding  Officer  and  Director 
Defense  Services  Administration 
Defense  Documentation  Center 
Cameron  Station,  Building  5 
5010  Duke  Street 
Alexandria,  VA  2231k 

Office  of  Naval  Research 
Resident  Representative 
Room  No.  582,  Federal  Building 
Austin,  TX  78701 

Electroacoustics  Group,  ARL/UT 

Electroacoustics  Division  I,  ARL/UT 

Computer  and  Marine  Sciences  Group,  ARL/lTT 

Signal  Physics  Division,  ARL/UT 

Sonar  Research  Group,  ARL/UT 

Sonar  Research  Division  HI,  ARL/UT 

Sonar  Engineering  Group,  ARL/UT 

Sonar  Engineering  Division  I,  ARL/UT 

Sonar  Engineering  Division  II,  ARL/UT 

Sonar  Development  Division,  ARL/UT 

Underwater  Systems  Division,  ARL/UT 


* 

\ 


maamm 


196 


p I ' 

Distribution  List  for  ARL-TR-76-2  (Revised)  (Cont'd) 


Copy  No. 

83 

Garland  R.  Barnard,  ARL/lTT 

6k 

Loyd  D.  Hampton,  ARL/lTT 

85 

Kent  M.  Mize,  ARL/lTT 

86 

Bruce  0.  Moses,  ARL/lTT 

87 

C.  Richard  Reeves,  ARL/lTT 

88 

Jack  A.  Shooter,  ARL/lTT 

89 

David  A.  Smith,  ARL/lTT 

90 

Reuben  H.  Wallace,  ARL/lTT 

91 

Library,  ARL/lTT 

92-100 

Reserve,  ARL/lTT 

